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ABSTRACT 

We calculate the coupled hydrostatic and ionization structures of spherically 
^ ; symmetric, pressure-supported gas clouds that are confined by gravitationally 

dominant dark-matter (DM) mini-halos and by an external bounding pressure 
. ^ ; provided by a hot medium. We focus on clouds that are photoionized and heated 

by the present-day background metagalactic field and determine the conditions 
for the formation of warm (WNM), and multi-phased (CNM/WNM) neutral 
atomic hydrogen (HI) cores in the DM-dominated clouds. 

We consider ACDM dark-matter halos with cuspy (NFW) and constant den- 
sity (Burkert) cores. We compute models for a wide range of halo masses, total 
cloud gas masses, and external bounding pressures. 

We present models for the pressure-supported HI structures observed in the 
Local Group dwarf irregular galaxies Leo A and Sag DIG. We find that the hydro- 
gen gas becomes neutral for projected HI column densities exceeding 10 19 cm -2 . 
We identify the HI cloud boundaries observed in Leo A and Sag DIG with the 
ionization fronts, and we derive an upper limit of Pmu/k ~ 100 cm~ 3 K for the 
ambient pressure of the intergalactic medium in the Local Group. The observed 
HI gas scale heights in Leo A and Sag DIG imply characteristic DM densities of 
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1.2 amu cm -3 (or 0.03M Q pc~ 3 ), consistent with the DM densities previously in- 
ferred via HI rotation curve studies of dwarf and low-surface brightness galaxies. 
Leo A and Sag DIG obey the scaling correlations that are expected for typical 
(median) DM halos in a ACDM cosmology, provided the halos contain constant 
density cores, as suggested by Burkert. 

We construct explicit "mini-halo" models for the multi-phased (and low- 
metallicity) compact high- velocity HI clouds (CHVCs). If the CHVC halos are 
drawn from the same family of halos that successfully reproduce the dwarf galaxy 
observations, then the CHVCs must be "circumgalactic" objects, with character- 
istic distances of 150 kpc. For such systems we find that multi-phased behavior 
occurs for peak WNM HI column densities between 2 x 10 19 and 1 x 10 20 cm -2 , 
consistent with observations. In contrast, if the CHVCs are "extragalactic" ob- 
jects with distances £ 750 kpc, then their associated halos must be very "un- 
derconcentrated" , with characteristic DM densities £j 0.08 cm -3 , much smaller 
than expected for their mass, and significantly smaller than observed in the dwarf 
galaxies. Furthermore, multi-phased cores then require higher shielding columns. 
We favor the circumgalactic hypothesis. If the large population of CHVCs repre- 
sent "missing low-mass DM satellites" of the Galaxy, then these HI clouds must 
be pressure-confined to keep the gas neutral. For an implied CHVC mini-halo 
scale velocity of v s = 12 km s -1 , the confining pressure must exceed ~ 50 cm -3 
K. A hot (~2x 10 6 K) Galactic corona could provide the required pressure at 
150 kpc. 

Our static mini-halo models are able to account for many properties of the 
CHVCs, including their observed peak HI columns, core sizes, and multi-phased 
behavior. However, important difficulties remain, including the presence in some 
objects of extended low column density HI wings, and Ha emission line fluxes in 
several CHVCs that are significantly larger than expected. 

Subject headings: Galaxy: general, formation, evolution - Galaxies: Local Group 
- ISM: high velocity clouds, intergalactic medium - Cosmology: dark matter 



1. Introduction 

High velocity clouds (HVCs) are atomic hydrogen (HI) clouds with radial velocities in- 
consistent with gas in differential circular rotation in the Galactic disk. Since their discovery 
via 21 cm observations by Muller, Oort & Raimond (1963) the HVCs have been observed, 
surveyed, and catalogued in ever increasing detail (Bajaja et al. 1987; Hulsbosch & Wakker 
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1988; Wakker k van Woerden 1991; Hartmann & Burton 1997; Braun & Burton 2000; Brims 
et al. 2000; Burton, Braun & Chengalur 2001; Putman et al. 2002; Lockman et al. 2002; de 
Heij, Braun & Burton 2002a). However, the nature and origin of the HVCs have remained 
the subject of considerable and unresolved debate (Oort 1966; Wakker & van Woerden 1997; 
Blitz et al. 1999; Wakker, van Woerden & Gibson 1999; Gibson et al. 2002; Sembach 2002). 

The HVCs consist of several hundred compact, kinematically distinct, clouds with typ- 
ical angular diameters of ~ 1 — 2° (Braun & Burton 2000; Putman et al. 2002), as well 
as several large (diffuse) and kinematically continuous complexes (e.g., the "A, B, C, H, 
& M" clouds) that extend over ~ 100 square degrees. The HVCs appear to be devoid of 
any stellar counterparts. The radial velocities of the HVCs identified in the Wakker & van 
Woerden (1991) catalogue range from —464 to +297 km s" 1 relative to the local standard 
of rest, de Heij, Braun & Burton (2002b) identified 67 distinct compact high velocity clouds 
(CHVCs) in the single-dish "Leiden/Dwingeloo" HI survey (Hartmann & Burton 1997). 
Putman et al. (2002) identified an additional 179 CHVCs in the (more sensitive) southern- 
hemisphere "HI Parkes All Sky Survey" (HIPASS). Most of the CHVCs are infalling. The 
"Leiden/Dwingeloo" survey, carried out with angular and velocity resolutions of 0?5 and 
1 km s -1 respectively, shows that the typical (one-dimensional) velocity dispersions of the 
CHVCs arc ~ 14 km s _1 with little scatter about this value (Blitz et al. 1999). The ve- 
locity dispersions imply gas temperatures ~ 10 4 K, consistent with a warm medium (WM) 
possibly consisting of a mixture of a warm neutral (WNM) and ionized (WIM) gas. Typi- 
cally, N H i ~5x 10 18 cm -2 averaged over the extents of the CHVCs. However, the CHVCs 
are barely resolved with the (25-meter) Dwingeloo and (64-meter) Parkes telescopes. More 
recently Burton et al. (2001) were able to resolve the spatial structures of ten CHVCs in 
higher resolution (3.5 x 3.7 arcmin) observations with the (305-meter) Arecibo telescope. 
They found that the typical observed (projected) 1/e exponential scale- lengths of the HI 
distributions equal 690 arcsec, with central HI columns ranging from 2 x 10 19 to 2 x 10 20 
cm -2 . In some objects, gas was detected out to ~ 1° from the cloud centers down to column 
densities of ~ 2 x 10 17 cm~ 2 . The line widths are slightly narrower in the Arecibo data set 
compared to the Leiden/Dwingeloo observations, with corresponding velocity dispersions of 
~ 11 km s -1 . The velocity gradients across the clouds are relatively small, ~ 10 km s -1 
degree -1 , and the WNM cloud dynamics are controlled primarily by thermal motions. 

In another important development Braun & Burton (2000) and de Heij, Braun & Burton 
(2002c) carried out high resolution Westerbork interferometric observations (with ~ 1 arcmin 
synthesized beams) and found that at least some CHVCs contain high column density (up to 
~ 10 21 cm -2 ) cores, with typical radii ~ 10 arcmin, and line profiles as narrow as 2 km s -1 . 
The narrow line widths imply HI gas in a cold neutral medium (CNM; T ;$ 100 K). Previous 
to this work CNM cores were known to exist only in the extended A, H, and M complexes 
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(Wakker & Schwarz 1991). The recent single-dish and interferometric observations reveal 
that "core/envelope" CNM/WNM structures may be a common feature of the CHVCs. 

The distance is the critical unknown for the HVCs. The large complexes are plausibly 
nearby objects, and stellar absorption-line observations in the large M and A clouds (Danly, 
Albert & Kuntz 1993; van Woerden et al. 1999) constrain their distances to within 5 and 10 
kpc respectively. However, no direct distance determinations are yet available for the large 
population of CHVCs. 

Many theories for the origin of the HVCs have been proposed over the years (Wakker & 
van Woerden 1997). Oort (e.g. 1966; 1970) recognized that if the HVCs are gravitationally 
bound objects with masses dominated mainly by the observed HI gas they must be distant, 
with distances d £ Av 2 f / nGm-RN-ftfl , where Av is the observed velocity dispersion, iVm is 
the mean HI column within angular radius 9, and / is the fraction of the total dynamical 
mass present as HI gas. If the clouds are bound, then the required value of / is proportional 
to the assumed distance d, since the HI mass varies as d 2 , whereas the dynamical mass 
M dyn = Av 2 d9/G oc d. For / = 1 the observations imply implausibly large "binding 
distances" d £ 10 Mpc, and HI masses M ffl = ird 2 9 2 m} i N iU >, 10 9 M Q . Oort's preferred 
interpretation was that the HVCs are transient "circumgalactic" ( £ 100 kpc) primordial gas 
clouds that have recently condensed out of the intergalactic medium, and are only now being 
accreted onto the disk in the end-stages of the Galaxy formation process. Shapiro & Field 
(1976) introduced the idea of the Galactic fountain (see also Bregman 1996) in which hot gas 
ejected from the disk cools and condenses into neutral clouds that fall back onto the Galaxy 
in a continuing cycle. Giovanelli (1981) suggested that the HVCs are scattered fragments of 
the Magellanic Stream. 1 A Local Group 2 (LG) origin for the HVCs was also proposed and 
discussed by various authors in the years immediately following their discovery (Oort 1966; 
Verschuur 1969; Einasto et al. 1974; Eichler 1976). In this picture the HVCs are genuine 
extragalactic objects associated with the LG system, and are at typical distances of ~ 1 



1 The "Magellanic Stream" (Mathewson, Cleary & Murray 1974) is an example of extended HVC gas 
that is relatively well understood. It consists of a narrow trail of HI gas that emanates from the Magellanic 
Clouds, and stretches over 100° across the sky. The stream is almost certainly gas that has been tidally 
(Gardiner & Noguchi 1996) or ram-pressure stripped (Moore & Davis 1994) from the Magellanic Clouds as 
these satellite galaxies orbit the Milky Way. 

2 We recall that the LG is a dynamically bound system of galaxies consisting of the Milky Way and its 
two largest neighbors M31 and M33, and 40 or more smaller and fainter dwarf galaxies (Mateo 1998). The 
LG is typical of the many small groups and clusters of galaxies situated along the outer boundary of the 
Virgo supercluster (van den Bergh 1999). The LG radius defined as the zero- velocity surface relative to the 
Hubble flow is ~ 1.5 Mpc, and the characteristic velocity dispersion is ~ 60 km s _1 (Sandage 1986). The 
mass of the LG is dominated by the Milky Way and M31, and equals about 3.3 x 10 12 M©. 
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Mpc, rather than being a structural feature of the Galaxy and its immediate surroundings. 

The distances to the clouds may be constrained via observations of optical Ha recom- 
bination line emission that provides a measure of the ionizing radiation field incident on the 
clouds (Weiner & Williams 1996; Bland-Hawthorn et al. 1998; Tufte et al. 1998; Weiner, Vo- 
gel & Williams 2001). At sufficiently large distances from the Galaxy, the dominant radiation 
field becomes the weak metagalactic background and the clouds should be weak Ha emitters. 
Nearby clouds might be ionized by more intense "leakage radiation" from the Galactic disk, 
and could be more intense Ha sources. Sensitive optical Fabry-Perot Ha observations have 
been reported and compiled by Weiner et al. (2001) for several HVC complexes and one 
isolated CHVC. Tufte et al. (2002) have reported Ha observations in five CHVCs. For the 
complexes Weiner et al. find that the measured Ha surface brightnesses range from 41 to 
1680 milli-Rayleighs 3 (mR), implying photoionizing Lyman continuum (Lyc) intensities of 
1.5xl0 4 to 6.0xl0 5 photons cm 2 s 1 sr 1 . Weiner et al. detected an Ha intensity of 48 
mR in CHVC 230+61+165 (an object included in the Burton et al. [2001] Arecibo data 
set). Tufte et al. (2002) report Ha intensities in the range of 20 to 140 mR in four CHVCs, 
and an upper limit of 20 mR in an additional CHVC. The detected fluxes are larger than 
the (2a) upper limit of 20 mR set by Vogel et al. (1995) for the Ha intensity of the isolated 
extragalactic Giovanelli & Haynes (1989) cloud that Vogel et al. argue is probing the meta- 
galactic field (see also Madsen et al. 2001; Weymann et al. 2001). Tufte et al. argue that the 
relatively large CHVC Ha fluxes imply that the CHVCs are located within the Galactic halo. 
However, the distance estimates are uncertain due to uncertainties in the actual strength 
of the metagalactic field, the Galactic radiation escape fraction, and the role of collisional 
ionization. 

The metallicity and dust content are additional potentially important clues to the origin 
of the HVCs. Thermal (~ 100 jum) emission signatures of dust grains have been searched 
for but not detected (Wakker & Boulanger 1986). However, a few metallicity measurements 
are available, mainly via ultraviolet absorption line studies toward background quasars and 
Seyfert galaxies. Significantly subsolar (though not purely primordial) metal abundances are 
indicated, primarily via observations of gas phase sulfur 4 . Low metallicities would appear 
to be inconsistent with a Galactic fountain, and favor models in which the gas originates 



3 1 Raylcigh = 10 6 /47r photons cm~ 2 s -1 sr -1 . For case-B recombination at a temperature of 10 4 K, a 
surface brightness of 1 Raylcigh is produced by a slab with emission measure EM = 2.8 cm~ 6 pc, which 
for two-sided illumination in ionization equilibrium, is produced by a Lyman continuum photon intensity 
J* = 3.6 x 10 5 photons cm~ 2 s _1 sr -1 . 

4 Sulfur is generally considered to be robust probe of the intrinsic metallicity since S remains undepleted 
in diffuse gas (Jenkins et al. 1987). 



- 6- 



outside the Galaxy. Lu et al. (1998) observed HVC 287+22+240 using the Goddard High 
Resolution Spectrograph (GHRS) on board Hubble Space Telescope (HST) and detected 
several SII and Fell absorption lines. Lu et al. derived S/H=0.25±0.07, Fe/H=0.033±0.006, 
and S/Fe=7.6 ± 2.2 (relative to solar abundances). They interpreted the supersolar S/Fe 
ratio as being due to significant iron depletion and therefore evidence for the presence of 
dust grains. In view of similar abundances observed in the Magellanic Clouds, as well as 
the proximity of 287+22+240 to these objects, Lu et al. concluded that this HVC is a 
"leading arm" component of the Magellanic stream. Wakker et al. (1999) carried out GHRS 
spectroscopy of HVC gas in complex C and inferred S/H=0.089± 0.024. They favor the Oort 
hypothesis, and interpret complex C as infalling primordial gas that has been "contaminated" 
by metals as it interacts with the Galactic halo. Braun & Burton (2000) derived a similar 
low metallicity for one of the CHVCs in their sample (125+41-207). They combined their 
HI observations together with Mgll absorption measurements reported by Bowen & Blades 
(1993), and inferred a Mg/H abundance in the range of 0.04 to 0.07 relative to solar. 

Recently, Murphy et al. (2000) and Sembach et al. (2000) carried out Far Ultraviolet 
Spectroscopic Explorer (FUSE) absorption-line studies of several HVCs. In complex C Mur- 
phy et al. found Fe/H ~ 0.5 (relative to solar) that may indicate undepleted iron and hence a 
very low dust content. Sembach et al. (2002) set a metallicity limit of O/H < 0.46 in CHVC 
224-83-197 via the detection of Ol absorption. An important result of the FUSE observa- 
tions was the detection of OVI absorption tracing highly ionized gas in several HVCs. These 
observations complement earlier GHRS observations of CIV, SilV and NV absorption in high 
velocity gas that Sembach et al. (1999) suggested is produced in the ionized envelopes of 
neutral hydrogen HVCs (as the ionized clouds do not appear to be 21 cm sources). Sembach 
et al. (1999) argued that the observed C +3 , Si +3 and N +4 abundances are consistent with 
low pressure (P/ k ~ 2 cm -3 K) gas photoionized by an assumed Lyc flux at a level plausibly 
consistent with the metagalactic background. They concluded that the CIV, SilV and NV 
absorbers are located at large distances and are Local Group objects. The OVI abundances 
are too large to be compatible with photoionization, and may be produced collisionally at 
the interfaces between the HVCs and a tenuous Galactic corona or Local Group medium 
through which the they are moving (Sembach et al. 2000). We note that collisional ioniza- 
tion (in gas with pressures considerably larger than 2 cm~ 3 K) may also contribute to the 
production of the C +3 , Si +3 and N+ 4 ions. 

Evidence of an ambient medium in which the CHVCs may be embedded could also 
provide clues to their origin. Briins, Kerp & Pagels (2001) mapped one of the CHVCs 
(125+41-207) in the Burton et al. sample with the 100m Effelsberg telescope. They found 
that the WNM component is asymmetric with a "cometary" appearance. The head-tail 
morphology suggests that gas is being stripped off of the main body of the cloud as it moves 
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through an enveloping medium. 

The Local Group hypothesis for the HVCs, particularly the CHVCs, has received re- 
newed interest recently for additional observational and theoretical reasons. Bajaja et al. 
(1987) noted that asymmetries in the "position- velocity" diagram of the HVC ensemble 
could be understood if the HVCs are moving within the LG as a whole. More recently 
Blitz et al. (1999) and Braun & Burton (1999) demonstrated that the dispersion in CHVC 
velocities is minimized if the cloud velocities are measured relative to the LG barycenter 
(see however, Gibson et al. 2000). Furthermore, Blitz et al. presented a simulation of the 
LG dynamics that shows that "test particles" interacting with the Milky Way, M31, and 
the distant Virgo cluster, are (during a Hubble time) drawn into an elongated filamentary 
structure that reproduces the observed distribution and velocities of HVCs in the directions 
of the LG barycenter and anti-barycenter (see however Moore & Putman 2001; de Heij et 
al. 2002a). These kinematic considerations are strong evidence for a LG origin, although 
the cloud distances are not very well constrained. Blitz et al. concluded that the typical 
distances of the CHVCs are d ^ 750 kpc. For d = 750 kpc, a characteristic radius of ~ 2.5 
kpc and a central gas density of ~ 6 x 10~ 3 cm~ 3 are implied by the Burton et al. (2001) 
observations. The HI mass in each cloud is ~ 1 x 1O 7 M , and the total HI mass in the 
system of CHVCs is £ 10 9 M Q . Such large masses may be in conflict with the apparent 
lack of similar objects around other galaxies and groups (Zwaan & Briggs 2000; Rosenberg 
& Schneider 2002). If the CHVCs are long-lived, gravitationally bound objects as conjec- 
tured by Blitz et al., then the measured velocity dispersions imply total dynamical masses 
M dyn k 5 x 1O 8 M , comparable to dwarf galaxy masses. In this picture the HI mass is a 
fraction / = 0.02 of the total dynamical mass, and the 3.4 kpc radius corresponds to the 
"scale height" ~ GMd yn /At> 2 of the gravitationally confined gas. Blitz et al. suggested that 
the dynamical mass is dominated by dark-matter and that the CHVCs are tracing individual 
dark-matter "halos". 

The notion that the CHVCs represent dark-matter halos is motivated by the cosmolog- 
ical theory of hierarchical structure formation. In this theory, bound systems ( "halos" ) of 
non-baryonic, dark matter evolve in a hierarchical collapse, initiated by small gravitation- 
ally unstable primordial density fluctuations (Gunn & Gott 1972; Press & Schechter 1974; 
Blumenthal et al. 1984; Navarro, Frenk & White 1997; Bullock et al. 2001). Low-mass halos 
virialize first, and then merge into progressively more massive systems. Galaxies form as gas 
accumulates, cools, and collapses inside the virialized DM halos (White & Rees 1978). De- 
tailed numerical simulations of the evolving collapse of dissipationless cold (non-relativistic) 
dark matter (CDM) have been carried out by several investigators (e.g. Navarro et al. 1997; 
hereafter NFW). The simulations predict the existence of many more low-mass halos than 
observed at the faint-end (i.e. dwarf galaxy scale) of the galaxy luminosity function. Various 
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suggestions have been made to account for this discrepancy, including: blow out of gas by 
winds and supernovae produced in small initial starbursts (Dekel & Silk 1986; Tegmark, Silk 
& Evrard 1993); photoionization heating of the baryonic component by the metagalactic 
radiation field, that inhibits the collapse of gas into the low-mass halos (Efstathiou 1992; 
Thoul & Weinberg 1996; Kepner, Babul & Spergel 1997; Kitayama & Ikeuchi 2000; Bullock 
Kravtsov & Weinberg 2000) and expulsion of gas trapped within the low-mass halos during 
the epoch of "reionization" (Barkana & Loeb 1999). In all of these scenarios the low-mass 
halos (or "mini-halos" as coined by Rees 1986) continue to exist, but because star-formation 
has been inhibited in them, they are not observable as galaxies. 

In two recent studies, Klypin et al. (1999) and Moore et al. (1999) presented high- 
resolution simulations of the evolution of dark matter halos at the mass and size scales of 
the Milky Way and M31 galaxies that dominate the LG. They found that many low-mass 
sub-halos survive and continue to accrete into each of the dominating galaxy halos to the 
present day. In particular, the predicted number of low-mass satellite halos appears to exceed 
the number of observed satellite dwarf galaxies within ~ 300 kpc of the Milky Way. Klypin 
et al. suggested that the CHVCs are in fact the "missing DM satellites" that appear in their 
simulations. In this picture the CHVCs are "circumgalactic" mini-halos associated with the 
Milky Way, as opposed to "extragalactic" objects associated with the LG as proposed by 
Blitz et al. The presence of DM substructure may be detectable via the gravitational lensing 
of background quasars (Mao & Schneider 1998; Metcalf & Madau 2001; Dalai & Kochanek 
2002). 

The CNM/WNM structures observed in the CHVCs are in many ways similar to the 
multi-phased HI distributions observed in gas-rich dwarf irregular (dlrrs) galaxies 5 . As 
summarized by Mateo et al. (1998), the observed HI masses in the dlrrs range from 0.1 
to 0.5 of the total galaxy mass (including the dark matter). The HI gas is generally much 
more extended than the "Holmberg radii" (0.1 — 1 kpc) of the stellar components as traced 
in the optical. Young & Lo (1996; 1997) carried out detailed Very Large Array (VLA) 
interferometric studies of the dlrrs Leo A and Sag DIG and found that the observed HI line 
profiles can be separated into broad and narrow components. The broad (WNM) components 
appear as "envelopes" throughout the mapped regions of the galaxies, whereas the narrow 
(CNM) components are associated with small clumps or cores located near star-forming 
regions. The WNM envelopes appear to be pressure supported as opposed to rotationally 
supported. In these respects the HI structures in the dwarfs are similar to those observed 
in the CHVCs. However, there are important differences (in addition to the fact that star 



5 Blitz & Robishaw (2000) have recently shown that some dwarf spheroidals (dSphs) may also be gas-rich, 
with the gas being present at large offsets from the optical stellar light. 
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forming regions are absent in the CHVCs). First, the peak HI column densities of the WNM 
gas in the dwarfs (up to a few 10 21 cm' 2 ) are significantly larger than in the CHVCs (up 
to ~ 10 20 cm -2 ). Second, the angular scale sizes of the HI distributions are smaller in the 
dwarfs, ~ 70 arcsec in Sag DIG, and ~ 160 in Leo A, compared to the CHVCs where the 
characteristic size is ~ 690 arcsec. These facts will be important in our analysis of the 
CHVCs and dwarf galaxies as possibly related objects. Another key difference is that the 
distances to the dwarfs are well determined (via stellar photometry) whereas the distances 
to the CHVCs are unknown. 

In this paper we quantitatively examine the hypothesis that the CHVCs are stable HI 
clouds confined by the gravitational potentials of dark-matter halos. We wish to determine 
whether viable "mini-halo" models for the CHVCs can be constructed as either "circum- 
galactic" (d ~ 100 kpc) or "extragalactic" (d ~ 1 Mpc) LG objects. Such clouds will be 
exposed to the background metagalactic radiation field that ionizes and heats the gas. The 
dark-matter confined clouds may also be subjected to bounding pressures, provided for ex- 
ample, by a hot circumgalactic Galactic corona, or by an intergalactic medium filling the 
Local Group. 

We address several questions. Given the observed HI column densities and distribu- 
tions what are the halo virial masses and characteristic dark matter densities required to 
confine the gas clouds? Are the required halo properties consistent with those expected from 
cosmological structure formation simulations? Under the assumption of photoionization by 
the metagalactic field, what are the minimum gas masses and column densities required for 
the formation of neutral hydrogen cores within the dark-matter halos? Are the observed HI 
distributions consistent with photoionization by the metagalactic field? Under what con- 
ditions do the cloud cores become multi-phased? Can halo models be constructed for the 
pressure supported HI clouds in the Local Group dwarf galaxies? If the CHVCs are the 
starless "cousins" of the dwarf galaxies, and the CHVC and dwarf galaxy halos are similar, 
what are the implied distances to the CHVCs? 

We construct general purpose models in which we compute the gas density distributions 
and ionization structures of pressure supported hydrostatic hydrogen gas clouds that are 
trapped in the potential wells of DM halos, are simultaneously heated and photoionized by 
an external radiation field, and are subjected to specified bounding pressures. Computations 
along these lines have been presented in the literature (with varying degrees of sophistication) 
mainly in the context of the mini-halo model for the intergalactic Lya clouds (e.g. Rees 
1986; Murakami & Ikeuchi 1990; Miralda-Escude & Rees 1993; Kepner et al. 1997). Of most 
relevance here are the computations presented by Kepner et al. (see also Corbelli & Salpeter 
1993) who studied the phase transitions (ionized/neutral/molecular) of hydrogen gas clouds 
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in DM halos. Their main focus was to show that low-mass halos form shielded neutral 
gas cores (a precondition for star-formation) only at late (z 1) cosmic times when the 
metagalactic photoionizing radiation field has become sufficiently diluted (Efstathiou 1992). 
However, Kepner et al. made several simplifying assumptions that make it difficult to apply 
their models to either the CHVCs or the LG dwarf galaxies. These included assuming a 
fixed relation between the halo mass and the characteristic DM density (see §2), adopting an 
arbitrary functional relationship between the total gas mass and halo mass, and the neglect 
of external pressure. In our analysis we relax all of these assumptions. 

In addition, we determine the conditions required for the formation of thermally unstable 
multi-phased cores. To do this we employ and extend the methods presented by Wolfire et 
al. (1995a) in their study of the multi-phased Galactic disk, and calculate the thermal 
equilibrium properties of low metallicity HI gas heated by the metagalactic field. We then 
incorporate the results of these computations into our dark-matter halo models. We note 
that Wolfire et al. (1995b) (see also Ferrara & Field 1994) carried out an explicit analysis of 
the HVC complexes observed by Wakker & Schwarz (1991). In particular, the cloud distances 
were constrained by arguing that the CNM detected by Wakker & Schwarz is possible only 
when the ambient gas pressure, exceeds a well defined minimum P m i n . In this analysis the 
gas pressure is provided by a Galactic X-ray corona. However, Wolfire et al. did not consider 
dark-matter dominated clouds. In such systems the central cloud pressures can become much 
larger than the ambient pressures. The transition from the WNM to the CNM phase then 
becomes possible in the cloud cores even when the ambient pressure is much smaller than 
p . 

1 mm- 

We first discuss in §2 the basic parameters and properties of the dark-matter halos we 
consider in our study, including also a discussion of the relevant cosmological scaling relations. 
In §3 we then describe our model clouds and discuss the methods we use to compute the cloud 
hydrostatics, radiative transfer, and HI phase structure. In §4 we discuss our representation 
of the metagalactic background field. In §5 we derive several useful analytic formulae for the 
gas distributions in DM halos. In §6 we present the results of our numerical computations. 
In §7 we model the HI gas structures in the dwarf galaxies Leo A and Sag DIG. In §8 
we consider models of the CHVCs as either circumgalactic or extragalactic mini-halos. A 
summary is presented in §9. 

2. Dark Matter Halos 

The analytic theory (Gunn & Gott 1972) of non-linear ("spherical top-hat") gravita- 
tional collapse shows that individual DM halos are most simply characterized by two inde- 
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pendent parameters (e.g. the characteristic density and halo radius). Numerical simulations 
of the dissipationless collpase of cold dark matter (NFW; Bullock et al. 2001) show that the 
DM distributions within the halos can be well represented by a spherically symmetric uni- 
versal density profile (the "NFW profile" ) over a wide range of mass scales. The simulations 
also show that for any particular cosmological model the two halo parameters are highly 
correlated. The correlation exists because a given mass-scale is associated with a particular 
virialization red-shift, with narrow dispersion. 

An important feature of the NFW profiles is that they are cuspy with diverging central 
DM densities. Central cusps in CDM halos are a predicted theoretical result of the dissi- 
pationless collapse. However, HI rotation curves in dark-matter dominated dwarfs and in 
several low surface brightness galaxies imply the presence of constant density cores (Moore 
1994; Burkert 1995; Blais-Ouellette, Amram & Carignan 2001; Firmani et al. 2001), although 
this conclusion may be hampered by beam smearing effects (van den Bosch & Swaters 2001; 
van den Bosch et al. 2000). In view of the HI rotation curves, Burkert (1995) introduced a 
semi-empirical profile (the "Burkert profile" ) in which the cuspy NFW core is replaced with 
a uniform density core. In our analysis we consider both NFW and Burkert profiles. 

Realistic dark matter halos may be complex triaxial structures, containing embedded 
clumps (sub-halos) and elongated filaments. In our analysis we represent the halos assuming 
they are spherically symmetric and smooth. We begin by defining the basic halo parameters 
for spherically symmetric systems, and we then discuss the cosmological correlation relation 
between the halo scale parameters. 



For spherically symmetric halos, the dark-matter density distribution may be written 

as 



where pd s is the characteristic "scale density" of the halo, f p is the density profile, and 
x = r/r s is the dimensionless radial coordinate, where r is the distance from the halo center 
and r s is the "scale radius" of the DM distribution. For realistic density profiles, the DM 
density is large for x < 1 and becomes small when x ^> 1. 

The enclosed DM mass at any radius is then 



2.1. 



Halo Parameters 



Pd = Pdsfp(x) 



(1) 



M d = M ds f M (x) 



(2) 
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where 

f M (x) = 3 [ X f p (x')x' 2 dx' , (3) 
Jo 

and 

47T 

M ds = yp rfs r s 3 , (4) 
is the characteristic or "scale mass" of the halo. 

The gravitational potential of the halo may be written as 

V? = v 2 Mx) (5) 

where 

, ( x f x fM(x') A , 

U( x ) = I —J^ dx > ( 6 ) 



x 

and where we define the "scale (circular) velocity" 

2 GM ds 4vrG 2 

^ s = = —z-Pdsr a , (7) 

r s o 

where G is the gravitational constant. 

The circular velocity at any distance from the halo center is 

v = v s f v (x) (8) 

where 

/„ W = [M£)1 I/2 . (9) 

L X J 

It is evident from equations (l)-(3) that for an assumed density profile f p the physical 
scales of a halo are fixed by two independent structural parameters. These may be chosen 
as two out of the four mutually dependent parameters pd s , r s , v s or M^ s . E.g., a halo may 
be defined by the scale density and radius, pd s and r s (using eq. [1]), or by the scale velocity 
and radius, v s and r s (using eq. [5]). The functions f p , Jm, and f v for NFW and Burkert 
halos are listed in Table 5 in Appendix A. 



2.2. Cosmological Relations 

In cosmological studies the two independent halo scale parameters are often given in 
terms of the virial radius, r vir (or the virial mass M vir ) , and the "concentration parameter," 
x vir = r vir /r s (often denoted by "C"). The virial radius is defined as the radius within which 
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the mean dark-matter density of the halo exceeds the background matter density of the 
universe, p u , by a factor A ~ 200, corresponding to the overdensities expected in collapsed 
and virialized objects (Gunn & Gott 1972). The virial mass, 

M vir = yAp u r 3 ir (10) 

is the DM mass enclosed within the virial radius. Thus, M vir = M ds f M (x v i r ). 
With these definitions it follows that 

/jW^vir) 



Pds = A Pu -^- , (11) 



1/3 m 1 

MlH— , (12) 



AnA Pu J vir x 



vir 



and 



4tt „ q . \ V6 



G 3 A Pu M l J:\— ^1 (13) 

In these expressions, the matter density p u and the overdensity factor A are fixed by the 
underlying cosmological model. The halo scale parameters p s , r s , and v s , are then determined 
by a choice of M vir and x vir . 

An important feature of the CDM collapse simulations is the finding that at any cosmic 
epoch the two independent halo scale parameters are correlated, over a wide range of halo 
masses (Navarro, Frenk & White 1996, 1997; Bullock et al. 2001). The correlation can be 
expressed in several ways, e.g., as a relation, x V i r (M V i r ), between the concentration parameter 
and the virial mass, or as a relation, r s (pd s ), between the scale radius and scale density. 
The correlation exists because a given mass scale (as defined for example by r s and pd s ) is 
associated with a narrow range of virialization redshifts, while the halo scale density pd s is 
effectively determined by the value of the cosmic matter density at the epoch of virialization. 
The numerical simulations show that for any assumed cosmological model the correlation 
is well approximated by a simple power-law relation x v - a ~ M" ir , where for low masses the 
index a is a weak function of the mass. For a given virial mass a dispersion of concentrations 
is found, with increasing departures from the typical concentration being less likely. 

We adopt the results of Bullock et al. (2001) for simulations carried out assuming the 
"concordance" ACDM cosmological model (Bahcall et al. 1999). In this model the matter 
and vacuum density parameters Q m = 0.3 and Q\ = 0.7, the Hubble constant H = 70 km 
s _1 Mpc" 1 , the root-mean-square (rms) amplitude of mass fluctuations in spheres of radius 
8/2T 1 Mpc er 8 = 1, and the present day matter density of the universe is p u = 2.76 x 10~ 30 



3 



1/2 
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g cm 3 = 1.66 x 10 6 amu cm 3 . 6 For this model the overdensity factor A for virialized 
halos equals 340 (Bryan & Norman 1997). The fraction of matter contained in baryons is 
Vl b /Q m = 0.088 (assuming Q b h 2 = 0.013). 

We are interested in the predicted halo properties at the current epoch (z — 0). Bullock 
et al. carried out numerical simulations for halos with present-day virial masses ranging from 
10 14 M down to 10 11 M Q . They also presented a "toy model" for the halo evolution that 
accurately reproduces the relation between concentration parameter and virial mass found 
in the collapse simulations. The numerical results show that the dispersion in concentrations 
(for a given halos mass) is well approximated by a log-normal distribution in which a standard 
deviations correspond to departures A(loga; v i r ) = 0.14a from the median concentration 7 . 
For halos between 10 11 and 10 14 M , Bullock et al. find that the expression 

^lasxicF^^p (14) 

provides an accurate fit to the results of their numerical simulation and associated toy model. 

In our study we are interested in low- mass halos with virial masses in the range ~ 10 s 
to 10 10 M , at and somewhat below the "dwarf galaxy" mass scale. This mass range is well 
below the resolution of the numerical simulations presented by Bullock et al. We therefore 
use their ACDM toy model (Bullock priv. comm.) to compute the concentrations for low- 
mass halos. In Figure 1 we plot the computed concentrations for median (a = 0) halos with 
virial masses equal to 10 8 , 10 9 , 10 10 and 10 n M o . We find that the fit 

/ M ■ \ -°- 08 

** = 27x 100 (to) (15) 

is accurate to within 3% of the toy-model computations for median halos. In this expression 
we have also assumed that, as for the high- mass halos, one-sigma deviations correspond 
to A(logrr vir ) = 0.14 departures from median halos. In Figure 2 we plot the x v i r (M vir ) 
relation, as given by equation (15) for families of +3<r overconcentrated, median, and — 3cr 
underconcentrated halos. The weak mass dependence, and the slight flattening of the power- 
index, from a = —0.13 in equation (14) to a = —0.08 in equation (15), is consistent with the 
behavior of the (post-recombination) cold dark matter fluctuation spectrum, which remains 
flat at low masses (Peebles 1982; Blumenthal et al. 1984; Eke, Navarro & Steinmetz 2001). 



6 A density of 1 amu cm" 3 equals 1.66 x 1CT 24 g cm" 3 = 0.93 GeV cm" 3 = 2.45 x 10 _2 M Q pc~ 3 . 

7 Due to a typographical error the la deviations are incorrectly given as A(logx v i r ) = 0.18 in Bullock et 
al. (2001) (Bullock priv. comm.). 
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For a given value of a, equation (15) together with equations (11), (12), and (13) may 
be used to compute all of the halo properties as functions of a single halo parameter. As 
is indicated by these equations, the resulting scale parameters depend on the enclosed mass 
function /m(^)- We find that virtually identical results are obtained for NFW halos and 
Burkert halos. In Figure 2 we plot the halo concentration, x v - ir , the virial mass, M v ; r , the 
scale density, rids = pds/n^u (amu cm -3 ), and the halo scale radius, r s (kpc), as functions of 
the halo and scale velocity, v s (km s" 1 ) 8 . We display these quantities for +3<r, median, and 
— 3a halos. For median halos, the scale radius r s ranges from 0.29 to 6.4 kpc, for v s between 
10 and 100 km s -1 . The scale density increases with decreasing v s , reflecting the earlier 
collapse epoch of lower mass halos. However, as expected, the density dependence is weak. 
For median halos, rids decreases from 2.7 to 0.55 amu cm -3 for v s between 10 and 100 km 
s _1 . For a narrow range of scale velocities, large changes in the characteristic DM densities 
are possible only by varying a, i.e., by moving to over-concentrated or underconcentrated 
halos. 

In his study of rotationally supported dwarf galaxies, Burkert (1995) found that the 
radii of the inferred (constant density) DM halo cores are correlated with the magnitudes of 
the observed circular velocities, as expected theoretically. In Figure 2 we plot the observed 
9 correlation between r s and v s for the four low-mass dwarf spiral galaxies (DDO 154, DDO 
170 NGC 3109, and DDO 105) in Burkert's sample (see his Fig. 2). It is evident that these 
objects lie precisely on the expected relation for median (a = 0) halos. The dwarf galaxy 
rotation curve data thus provide empirical support for the theoretically derived correlation 
relation given by equation (15). 

Convenient analytic scaling relations for the various halo parameters as functions of v s , 
may be developed by setting fu(x) ~ lAOx 1 / 2 in equations (11) and (13). This approxima- 
tion for f M is accurate to within 12%, for both Burkert and NFW halos, for 4.5 < x < 50. 
It then follows that 

x vn = 32.1 x 10°- 14 %; 6 a225 , (16) 
M vir = 5.55 x 10 7 x 10- ail2 X 3 6 191 M , (17) 
n ds = 2.71 x 10°- 373 ^; 6 °- 638 amu cm" 3 , (18) 



8 In cosmological studies, halos are often parameterized by the circular velocity at the virial radius, 
iVir = v s [/m (^virV^vir] 1 / 2 , or the maximum circular velocity w max = 0.8v s (for NFW and Burkert potentials). 
We prefer to use v s as it directly fixes the potential- well depth (see eq. [5]). 

9 Burkcrt (1995) actually plotted the correlation between r s and the observed circular velocity at r s . We 
have transformed the observed velocities to scale velocities using equation (8) and Burkert's profile. This 
gives v(r s ) = 0.618w s . 
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and 

r s = 0.29 x io- - 186 ^ 319 kpc , (19) 

where v s6 = v s /10 km s _1 . These expressions are plotted as the dashed lines in Figure 2. 
They provide excellent analytic representations (to within 1% accuracy) of the halo scalings 
for median and underconcentrated halos, and are less accurate (~ 30%) for over-concentrated 
halos. 

In our computations we adopt the rr vir (M vir ) correlation relation as given by expres- 
sion (15), and the resulting halo parameters displayed in Figures 1 and 2, as theoretically 
motivated guides to the properties of "median" low-mass halos. In our study, we also explic- 
itly consider halos that are substantially "over-concentrated" or "under-concentrated" with 
respect to the median halos. 



3. Clouds 

We are interested in determining the structural properties of hydrostatic gas clouds 
that are confined by gravitationally dominant dark-matter halos, and are in ionization and 
thermal equilibrium with an external source of Lyman continuum (Lyc) and far-ultraviolet 
(FUV) radiation. We focus on ionization and heating by the metagalactic field. We envision 
clouds consisting of a "warm medium" (WM) with up to three components; (a) an outer 
layer of warm photoionized gas (WIM), that serves as a shielding layer for (b) a neutral core 
of warm HI gas (WNM), and possibly (c) a still smaller region within which the HI can exist 
as a multi-phased (WNM/CNM) mixture of warm and cold gas. When ionizing radiation is 
present, the WIM component will always be present. However, the WNM and WNM/CNM 
components will not always exist, and one of our goals is to determine the specific conditions 
required for their formation. 

We consider the possibility that the halo-cloud systems are embedded within a low- 
density hot ionized medium (HIM), provided for example by a Galactic X-ray corona in the 
circumgalactic environment, or an intergalactic medium in the Local Group environment. 
The HIM exerts a bounding pressure, Phim, on the outer surface of the WM cloud and serves 
to further confine the WM gas, particularly in the limit of weak DM halo potentials. We 
assume that the HIM is sufficiently hot such that it is completely unbound to any of the halos 
we consider. We also assume that the HIM is optically thin to the background radiation 
field. 

We wish to address several questions. For a given dark-matter halo and a total mass 
of warm gas, what are the radial extents and associated column densities of the HI gas 
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distributions? For small gas masses we expect the WM clouds to be photoionized throughout. 
In such WIM dominated clouds, any HI is present as a small neutral fraction within the 
ionized gas. As the gas mass is increased the WIM component can be expected to become 
optically thick to the ionizing radiation, leading to the formation of a neutral hydrogen 
WNM core. For an assumed halo and background field, what is the minimum WM mass 
required for the formation of a neutral WNM core? When WNM cores exist, how large 
are they? How do the core sizes and associated HI column densities depend on the halo 
scale parameters? If the hydrostatic pressure gradients are sufficiently large the central gas 
pressures may become large enough to allow a multi-phased WNM/CNM medium (or a 
complete conversion to CNM) in the neutral cores. For clouds in which the ionization and 
heating is dominated by the metagalactic field, when is multi-phased behavior possible in 
the neutral cores? 

To address these questions, we first construct models in which we solve the coupled 
equations of hydrostatic equilibrium and radiative transfer for warm spherically symmetric 
clouds that are (a) confined by gravitationally dominant halos, (b) subjected to a bounding 
pressure, and (c) exposed to an external source of Lyc and FUV radiation that ionizes the 
gas. We focus on irradiation by the metagalactic field. Because the temperatures of the warm 
components are expected to lie in a narrow range (~ 10 4 K) we specify the gas temperatures 
for the WIM and WNM components in advance in these computations. Given the resulting 
radial profiles for the hydrostatic gas pressures, densities, and fractional ionizations, we 
compute various quantities of observational interest including the HI gas masses and column 
densities, gas scale heights, locations of the ionization fronts (if present), cloud emission 
measures, Ha surface brightnesses, and other quantities of interest. Our treatment of the 
cloud hydrostatics and radiative transfer is described in §3.1. 

We then determine the phase states of the HI gas in the neutral cores (when such cores 
exist). The thermal phase state of the neutral gas at each cloud radius is determined by 
comparing the local value of the hydrostatic pressure P(x) to the local values of two critical 
pressures P m i n (x) and P m ax(x). At low pressures P < P m i n the gas is WNM, at intermediate 
pressures P mia < P < P max the gas is multi-phased, and at high pressures P > P max the gas 
must convert to CNM. 10 . 



10 The fact that interstellar HI gas can exist as WNM and CNM has been the subject of many theoretical 
investigations (e.g., Field, Goldsmith & Habing 1969; Zeldovich & Pilkelner 1969; McKee & Ostriker 1977; de 
Jong 1977; Drainc 1978; Ferriere, Zweibcl & Shull 1988; see also the reviews by Field 1975, and Kulkarni & 
Heiles 1987, 1988). Recent comprehensive theoretical studies have been presented by Wolfire et al. (1995a; 
2002). The two distinct phases are possible because the gas heating rates are generally proportional to 
the gas density, tin, whereas the gas cooling rates are proportional to n^. Thus, at low densities cooling 
becomes inefficient, and the temperature is driven to maximal (WNM) values of ~ 10 4 K set by Lya cooling. 
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As we discuss in §3.2 for a given incident radiation field the critical pressures P m i n and 
Pmax for shielded HI gas can be expressed as functions of the primary ionization rate ( p , 
independent of the cloud geometry (e.g. spherical vs. plane parallel). In §3.2 we present the 
results of thermal equilibrium computations which yield the required -P m in(Cp) and P ma x(Cp)- 

Given P mm (C P ) and P max (Cp) we determine the phase states of the HI gas in our halo- 
cloud models. In these models the local hydrostatic pressure depends on the depth of the 
halo potential well and the value of the bounding pressure, and increases from the surface to 
the core. Conversely, P m i n and P max decrease from the surface to the core as the ionization 
rate ( p decreases with increasing cloud depth and opacity. Thus, any neutral gas may be 
expected to be pure WNM at large radii where P < P mm , but may become multi-phased or 
fully converted to CNM at small radii where P min < P < P max or P > P max . For sufficiently 
deep halo potentials the conversion to CNM may occur even if the ambient bounding pressure 
is negligibly small. 

We determine the conditions for which P m i n (0) < P < P ma x(0) where P is the central 
pressure. In such clouds the multi-phased cores will extend out to a radius x where P(x) = 
Pmin(x). We assume that in the multiphased zones the CNM occupies negligible volume 
but remains uniformly mixed with the WNM. When P > P max (0) the HI in the cores 
will be driven entirely into the CNM phase. When this happens, thermal support is lost, 
and our hydrostatic solution in which warm gas is assumed to be present from the cloud 
center to the surface, becomes inconsistent with gas heating by the external field. Either 
additional (internal) heat sources must be postulated, or the cloud must collapse. Thus, 
in our computations we find the maximum masses, MwMmax, of warm gas that can be 
maintained by the external field within a given halo and an assumed bounding pressure. 
The total mass of CNM is unconstrained in our models, except that we assume it cannot 
dominate the gravitational potential. 



At high densities, cooling by fine structure line emission, primarily 158yum [CII] emission, becomes much 
more effective and the gas can cool to CNM temperatures ;$ 100 K. At intermediate densities a thermal 
instability occurs as the dominant cooling mechanism switches from Lya to [CII] emission. The instability 
and multiphased behavior occurs within a range of pressures - P m i n to P max - that depend on the detailed 
heating and cooling rates. For pressures less than P m ; n only the warm phase is possible, and for pressures 
greater than P max only the cold phase is possible. 
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3.1. Hydrostatics and Radiative Transfer 

In hydrostatic equilibrium 

dP = -pdip , (20) 

where P is the gas pressure, p is the gas density, and p> is the gravitational potential. The 
pressure P = pc 2 , where c 2 = kT / (m) is the isothermal sound speed of the gas, T is the gas 
temperature, and (m) is the mean mass per gas particle. We assume that the potential is 
dominated by the dark-matter halo so that (p = v 2 f ip (x), where v s is the halo scale velocity 
and f v is the dimensionless potential for either NFW or Burkert halos. It follows that the 
gas pressure P(x) is given by 

P(x) = Poexp [- f ^ll^ldx'} (21) 
L J c 2 g x' 2 J 

where /m(^) is the (dimensionless) DM mass distribution function, and Pq is the central 
gas pressure. It is evident that the hydrostatic structure is controlled by the ratio, v s /c g , of 
the halo scale velocity to the gas sound speed. The hydrostatic and ionization structures of 
the cloud are coupled since the sound speed c g depends on T and (m), which vary with the 
ionization state. 

In ionization equilibrium 

Cn(H°) = a B n(R + )n(e), (22) 

where ( is the total local photoionization rate, and n(H°), n(H + ) and n(e) are the atomic 
hydrogen, proton and electron densities respectively. We assume a case-B recombination 
coefficient «g. The local primary photoionization rate is 

POO 

Cp(x)=4tt/ a v J*(v,x)dv (23) 

J VQ 

where z/ is the Lyman limit frequency, a v is the photoabsorption cross section, and J*(u,x) 
is the local mean photon intensity 

l r 1 r -, 

J*(u,x) = -Jq{v) J exp\-a u N m (x,p)^dp, (24) 

where Jq(^) is the unattenuated isotropic radiation field (in units of photons cm~ 2 s _1 Hz -1 
sr _1 ), Nm(x,p) is the HI absorbing column from the cloud surface to x along a ray that is 
inclined by an angle 9 relative to the outward radial direction, and where p = cos#. 

For the range of hydrogen gas columns (up to ~ 10 22 cm' 2 ) that we consider in our 
computations, the gas opacity is dominated by a combination of hydrogen and helium pho- 
toabsorptions. We assume a primoridal helium abundance nu e /nn = 1/12 (Olive & Steigman 
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1995; Ballantyne, Ferland & Martin 2000) where nn and nu e are the total hydrogen and he- 
lium particle densities respectively. Photoabsorptions by heavier "metals" may be neglected 
for the low metallicities (Z £ 0.1) we assume (Morrison & McCammon 1983; Balucinska- 
Church & McCammon 1992). In our treatment of the radiative transfer we make the sim- 
plifying approximation that the hydrogen and helium ionization fractions are equal at all 
locations. We also assume that every helium recombination photon is absorbed in a hydrogen 
ionization. With these assumptions a v = a u ^+(l/12)a U; iie- We adopt the Balucinska-Church 
& McCammon expressions for ov,h an d o^Hc- In computing the total ionzation rate we also 
include the contributions of the secondary collisional ionizations produced by the primary 
photoelectrons (Shull & van Steenberg 1985). 

We construct models assuming that the background metagalactic field is the dominant 
source of radiation impinging on the clouds. Our representation of this field is displayed in 
Figure 3 in a v,J v plot. We discuss our fit to the metagalactic field in Appendix B. For our 
assumed field AnJ* = 1.28 x 10 4 photons s _1 cm~ 2 , where the Lyc integral 

AttJ* = J dVlJ^ du^- . (25) 

Most (93%) of the photons are emitted between 1 and 4 ryd. The unattenuated primary 
hydrogen ionization rate is C P ,o = 4 x 10 -14 s _1 . 

For our assumed helium abundance, 

P = ^(1 + x p )n u kT (26) 

where x p = n(H + )/nu is the proton fraction. The electron fraction x e = n(e)/nn = ff^V 
The mean mass per gas particle is 

/ v 16 1 

< m > = - 3 — p m » (27) 

where mn is the atomic hydrogen mass. 

In our computations we assume T = 10 4 K for both the WNM and the WIM. To verify 
this assumption for the WIM we have computed the gas temperature of low-metallicity pho- 
toionized gas using CLOUDY (version C90.04, Ferland et al. 1998) 11 assuming an ionizing 
field as given by our representation of the metagalactic field discussed in Appendix B. We 
find that the gas temperature of the ionized gas ranges from 8.0 x 10 3 to 1.2 x 10 4 K for gas 



11 Seewww.pa.uky.edu/~gary/cloudy/ 
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densities in the range 1CT 3 to 0.1 cm" 3 . The HI thermal phase computations we discuss in 
§3.2 show that the WNM gas temperature is close to 10 4 K. For 10 4 K gas, the associated 
gas sound speeds are c 9i wnm = 8.19 km s _1 for the WNM, and c 9i wim = H-6 km s -1 for the 
WIM. 

The local gas pressure, fractional ionization, and gas temperature, determine the local 
density of hydrogen nuclei riu via equation (26). We define the hydrogen gas (ionized plus 
neutral) density distribution / gas such that 

nu = n Hfi f gas (x) (28) 

where nn,o is the central density. We define the outer edge of the WM cloud as the radius 
xw/h at which the cloud pressure is equal to the pressure Phim of the enveloping hot ionized 
medium. Thus, £w/h is located at the interface between the WIM and the HIM. The total 
(hydrogen) cloud mass is then 

r x w/H 

M WM = 47rr 3 m H n H ,o / f gas {.x)x 2 dx . (29) 
Jo 

For a given halo potential and background field, a cloud model is selected by choosing a 
value for the total mass M W m given an assumed value for Prim, or equivalently, setting the 
boundary condition that the pressure at a selected value of xw/h equals Phim- 

We have constructed a code to solve the coupled hydrostatic and ionization equilibrium 
equations. In our numerical solution we specify Phim and £w/h, and then adopt an itera- 
tive procedure. First, we assume an ionization and opacity structure for the gas (initially 
by setting x e — 1 everywhere, and assuming that the cloud is optically thin). We then 
integrate equation (21) subject to the boundary condition P(xw/h) = Phim- This yields a 
hydrostatic pressure distribution P(x), and a value for the central pressure. P(x) and the 
assumed ionization structure determine the local values for the hydrogen gas densities. We 
then compute the local photoionization rates given the absorbing columns and opacities to 
the cloud surface. Given the ionization rates and the hydrogen gas densities we solve the 
equation of ionization equilibrium and find new values for the ionization fractions x e and 
cloud opacities. The local gas sound speeds are then updated, and the hydrostatic structure 
is recomputed etc. We iterate until the local ionization fractions have converged to within 
1%. 

Given the solutions to the hydrostatic and ionization structures, several quantities of 
interest may be computed, including the total gas column densities, the HI masses and 
columns, the cloud emission measures and Ha surface brightnesses, the gas scale heights, 
and the location of the ionization fronts. 
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The total column density of hydrogen nuclei along a line-of-sight with impact parameter 
x from the cloud center is 

N WM (x) = 2n R>0 r s / f gas (x') /2 _ 2]1/2 dx' , (30) 
and the HI column density along a line-of-sight a distance x from the cloud center is 
N m (x) = 2n R)0 r s / (1 - x p )f gas (x')————dx' , 

J x [X X \ I 

where O&p IS the local proton fraction. 

The emission measure along a line-of-sight with impact parameter x is 

EM(x) = 2n 2 lQ r s J x e x p flJx') _ /2 dx' (32) 

where x e and x p are the local electron and proton fractions. The associated Ha surface 
brightness in milli-Rayleighs may be written as 

J* a ( x ) = 10 9 'a eS EM = 364£M pc mR, (33) 

where a cS = 8.7 x 1CT 14 cm 3 s" 1 is the Ha effective recombination coefficient (Osterbrock 
1989) and EM pc is the emission measure in cm~ 6 pc. 

While formal solutions may be found for all values of Phim and xhim, several additional 
conditions must be satisfied for physically realizeable solutions. First, while the halos for- 
mally extend to infinite radius, we assume that effective physical boundaries exist beyond 
which both the baryons and dark-matter particles are subjected to external forces that read- 
ily detach them from the halo systems. As a plausible limit we assume that the radii of the 
halos and trapped clouds cannot exceed the virial radius. Thus, we require that x-w/h < ^vir 
(where x v ; r is the halo concentration parameter). If the halos are tidally truncated within 
the virial radius, then the limit becomes xw/h < r t/ r s, where r t is the tidal radius. 

Second, self-consistent solutions require that the DM mass dominates the gas mass at all 
radii, otherwise the assumption that the halo gravity dominates is violated. Self-consistent 
solutions must necessarily satisfy the condition M W m < ^dm(%/h)j where M dm (xw/h) is 
the dark-matter mass enclosed within the entire WM cloud. 

Third, the WM gas must be well bound to the halos. Otherwise, mass loss (especially 
from the outer WIM layers) will occur. We define bound WM clouds as those that satisfy 
the condition 

W _ v 2 s [f^(x w/H ) -U(oo)] 



1 2 C WIM 



> 1 (34) 
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where W is the gravitational potential energy per particle at xw/h, and T = (3/2)c^ IM is 
the mean thermal energy per particle in the outer WIM layer 12 . 

We conclude this section by making an important distinction between pressure- confined 
and gravitationally- confined clouds. Because the depths of the halo potential wells we are 
considering are finite (as x — > oo, j 9 — > 3 for NFW halos, and f v — > 3tt/A for Burkert halos, 
see Appendix A) a maximum possible pressure (density) contrast between large and small 
radii exists for a given value of v s /c g . When v s /c g is small the maximum pressure (density) 
contrast is small, and the gas extends to large halo radii x 1. In this limit a cloud of 
finite mass must be pressure-confined, and the halo gravity provides just a small pressure 
enhancement towards the cloud center. In pressure-confined clouds, properties such as the 
total cloud mass, or the internal gas pressure, are determined by conditions at the cloud 
surface via by the bounding pressure Phim- Pressure-confined clouds will tend to become 
unbound as the confining pressure becomes small. 

When v s /c g is large the pressure contrast becomes large, and the gas is effectively re- 
stricted to small halo radii i C 1. In this limit, which we refer to as the "small-x" limit, the 
cloud is gravitationally confined, and for any plausible value of the central pressure (or den- 
sity), the total gas mass remains finite for arbitrarily small bounding pressures. Clouds may 
be assumed to be in the small-x limit when v s /c g ^ 1.5 (see Appendix A). In gravitationally 
confined clouds, properties such as the total mass or internal gas pressure are determined by 
the central gas pressure, that is, by conditions at the cloud centers. Gravitationally confined 
clouds become increasingly bound as v s /c g becomes large and the potential wells deepen. 
When the observable gas is in the small-x limit, the implied dynamical mass Md yI1 = c g r gSuS /G 
(where r gas is the observed gas scale height) is a small fraction of the halo scale mass, and 
an even smaller fraction of the virial mass. In gravitationally confined clouds, much more 
dark-matter must exist beyond the observed gas structures, unless the halos are truncated 
at radii much less than their virial radii. 

In Appendix A we present simple analytic expresssions for the gas scale heights, gas 
masses, and column densities, for clouds in the small-x limit, and we derive an analytic 
condition for the transition from fully ionized clouds to clouds containg neutral cores. These 
analytic results (for both NFW and Burkert halos) are summarized in Table 6, and may be 
used as guides to understanding the numerical model results we present in §4, §5, and §6. 



12 A similar condition was adopted by Barkana & Loeb (1998) in their study of the photoevaporation of 
gas in low-mass halos at the reionization epoch. 



3.2. -Pmin &nd -P max 



We employ the methods described in Wolfire et al. (1995a, 2002) to compute P vs. nn 
"thermal phase diagrams" for low-metallicity HI gas heated by metagalactic radiation for 
a wide range of assumed HI shielding columns. We use the results of these computations 
to determine the critical pressures P min and P max as functions of the shielding column and 
associated primary ionization rate ( p . 

In our phase diagram computations the equations of thermal balance and of ionization 
and chemical equilibrium are solved self-consistently. In these computations we adopt a 
plane parallel geometry in treating the radiative transfer through the absorbing columns, 
and we assume an incident flux equal to 2irJ* where J„ is the mean intensity assumed in 
our spherically symmetric halo cloud models 13 . We calculate the hydrogen ionization rates 
using the same approximations for the Lyc absorption and opacity that we use in computing 
the ionization structure of the spherical halo clouds. 

The three major heating processes included in our HI phase computations are UV and X- 
ray ionization, FUV photoelectric emission, and molecular hydrogen (H 2 ) photodissociation. 
Their relative importance depends on the Lyc vs. FUV intensities, on the dust-to-gas mass 
ratio D/G, and the metallicity Z. The metallicity also controls the gas cooling efficiencies. 
In view of the generally low (but non-negligible) metal abundances observed in the HVCs and 
CHVCs we assume Z = 0.1, where Z = 1 corresponds to the Galactic (gas-phase) abundances 
of n /nu = 3.2 x 10" 4 for oxygen (Meyer, Jura & Cardelli 1998) and ric/nn = 1.4 x 10~ 4 
for carbon (Cardelli et al. 1996). The dust content and grain properties of the CHVCs 
are unknown, and we assume that the dust-to-gas ratio is proportional to the metallicity, 
and set D/G = 0.1. Following Wolfire et al. (2002) we adopt a slightly modified form of 
the Bakes & Tielens (1994) photoelectric heating efficiencies. The FUV field, that leads to 
photoelectric emission, produces the C + cooling ions and photodissociates H 2 , is attenuated 
by dust absorption. For D/G = 0.1 a photoelectric heating FUV optical depth of 1 occurs 
at an HI column density of 1.1 x 10 22 cm~ 2 (Roberge, Dalgarno & Flannery 1981). We 
assume that that H 2 forms on grain surfaces with a rate coefficient R = 3 x 10~ l7 Z cm 3 
s _1 proportional to the metallicity, and also via the standard gas phase sequence H+e~ — > 
H + v, H~+H + — >H 2 + v. The H abundance is limited by photodetachment, mainly 
by far-red photons, with a cross section that peaks near 8500 A (Wishart 1979). Detailed 



13 In spherical clouds the ionization rate near the surface of optically thick clouds is proportional to 2nJ* , 
whereas near the cloud center the ionization rate is proportional to At: J* (sec cq. [24]). In our plane-parallel 
computations we adopt an incident flux equal to a value of 2ttJ*, with the recognition that the resulting 
-Pmin and P max as functions of the ionization rate are insensitive to the assumed incident flux. 
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expressions for the various heating and cooling rates, and the basic atomic and molecular 
data, are given in Wolfire et al. (1995a) with some updates in Kaufman et al. (1999) and 
Wolfire et al. (2002). 

As a specific example we show in Figure 4 our numerical results for a shielding column 
Nni = 5 x 10 19 cm~ 2 . This particular shielding column is important since as we will show this 
is the typical peak HI column density for WM clouds in which the gas becomes multi-phased 
in the cloud cores. For our assumed representation of the metagalactic field (see Appendix 
B), and a shielding column of 5 x 10 19 cm -2 , the primary hydrogen ionization rate is equal 
to 1.1 x 10~ 18 s _1 . 

The resulting P vs. phase diagram is shown in Figure 4a for densities ranging from 
10~ 3 to 100 cm -3 . The thermally unstable region, defined as the zone where d(logP) / d(lognu) < 
0, is visible between the WNM and CNM branches, and a two-phased medium is possible 
for pressures between P m i n = 45.0 cm~ 3 K, and P max = 310 cm -3 K. The dominant heating 
and cooling rates per hydrogen nucleus, T and urA, are displayed in Figure 4b. UV/X-ray 
ionization dominates the heating of the WNM. The UV/X-ray heating rate decreases from 
2 x 10~ 28 to 2 x 10~ 29 ergs s _1 H _1 in the density range 10~ 3 to 100 cm~ 3 . The decline is 
due to a decrease in the fractional ionization. As x e decreases an increasing percentage of 
the UV/X-ray photoelectron energy is lost in atomic excitations and secondary ionizations 
rather than being deposited as heat. In contrast, the photoelectric heating rate (and effi- 
ciency) remains approximately constant, and is equal to ~ 2 x 10~ 29 ergs s _1 H _1 , across the 
entire density range. Consequently, the CNM is heated by a combination of UV/X-ray and 
photoelectric heating. At the highest densities the molecular formation rate becomes suffi- 
ciently large such that H 2 photodissociation heating dominates. The WNM is cooled by Lya 
emission, and the CNM is cooled by [CII] 158/xm fine structure emission. At high densities 
carbon remains neutral and [CI] 609/zm fine-structure line emission dominates the cooling. 
The computed fractional ionizations and gas temperatures are shown in panels Figures 4c 
and 4d. In the WNM, T ~ 10 4 K, and in the CNM, T = 10 — 100 K. 

In Figure 5a we display the P vs. n H phase diagrams for shielding columns ranging from 
1.0 x 10 18 cm~ 2 to 1.0 x 10 21 cm' 2 . As the shielding column increases, the ionization and 
heating rates decrease, and the critical pressures P max and P m in both decrease. In Figure 5b 
we plot P m in and P max as functions of the primary ionization rate, and the associated values 
of the shielding column. For our assumed field the ionization rate ranges from 8.9 x 10~ 16 
s _1 to 9.1 x 10~ 21 s _1 for shielding columns between 1.0 x 10 18 cm~ 2 and 3.0 x 10 21 cm~ 2 . 
For this range, P mm //c varies from 6.4 x 10 3 to 5.7 cm~ 3 K, and P max /k varies from 1.2 x 10 4 
to 11.2 cm" 3 K. 

For a fixed heating rate, P m i n and P max scale approximately as 1/Z (where Z is the 
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met alii city). This is because the fine-structure and metastable metal line cooling rates are 
proportional to Z, and the metal line cooling dominates at both P min and P max (see Fig. 4b). 

We use the computations of P m in(Cp) an d -Pmax(Cp) displayed in Figure 5b to determine 
the HI phase states of the gas within the halo clouds. As we have discussed, we compare the 
hydrostatic halo cloud pressure P to the values of P m i n and P max both of which vary with 
the local primary ionization rate within the halo clouds. When P < P m i n the gas is WNM, 
when P min < P < P max the gas is multi-phased, and when P > P max the gas must convert 
to CNM. 



4. Numerical Computations: M gas vs. M vir Diagrams 

We have carried out computations of the hydrostatic gas distributions, ionization struc- 
tures, and HI thermal phase properties, of dark-matter dominated clouds for a wide range of 
assumed total gas masses and halo masses. We consider both NFW and Burkert potentials, 
and a range of HIM bounding pressures. The computational results may be used to select 
models that fit the observed HI distributions and thermal phase properties observed in the 
Local Group dwarf galaxies and in the CHVCs. 

In Figures 6 through 11 we display a variety of computed cloud properties as functions 
of the total gas mass, Mwm, and halo virial mass M v ; r . We emphasize that the total gas mass 
Mwm refers to "warm gas" only, and includes the WIM and WNM, but excludes CNM. As 
discussed in §3, in computing the hydrostatic and ionization structures we adopt a metallicity 
of 0.1 Z Q , and we assume that T = 10 4 K in both the WNM and the WIM. We assume 
irradiation by the metagalactic field. We focus on the behavior of median (a = 0) Burkert 
halos and — 4<r NFW halos, since as we discuss in §7, for such halos the computed gas scale 
heights and WNM cloud sizes correspond to those observed in Local Group dwarf galaxies. 
For each family of halos we display results for HIM bounding pressures, equal to 1, 10, and 
100 cm" 3 K. 

The quantities shown in the M gas vs. M vir plots are; (a) the total cloud radius rw/H 
(kpc); (b) the total gas to DM mass ratio Mwm/Mdm within t w /h; (c) the radius twnm 
(kpc) of the WNM core (if present); (d) the total HI gas mass Mm (M Q ) contained in the 
WNM and WIM components; (e) the peak HI column density Nm (cm -2 ) integrated across 
the entire cloud diameter, and (f) the central gas pressure Po/k (cm~ 3 K). 

We recall (see §2) that for a given DM density profile (NFW or Burkert) and an assumed 
correlation relation, x vir (M vir ), between the halo concentration x vir = r vir /r s and the virial 
mass, a choice for M vir uniquely determines the halo structure and scale parameters. We 
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consider halos with M vir ranging from 10 s to 10 10 M . For median Burkert halos this mass 
range corresponds to scale velocities from 12 to 50 km s^ 1 . For the — 4<r NFW halos this 
mass range corresponds to v s ranging from 8 to 40 km s _1 . 

In Figures 6 through 11 the WM gas masses for which results are shown range from 
5 x 10 6 to 5 x 1O 7 M . This range was selected such that peak HI columns of 10 19 to 10 21 
cm" 2 , representative of the columns observed in the CHVCs and dwarfs, appear within the 
displayed M W m vs. M vir plots. 

For the purposes of the discussion that follows we refer to halos approaching ~ 10 8 M Q 
as "low-mass halos", and halos approaching ~ 1O 1O M as "high-mass halos". 

The total cloud radius r w /H is located at the outer edge of the WIM envelope where the 
cloud pressure equals the bounding pressure Phim- For fixed Phim, and a given halo mass 
M vir , the WM cloud radius increases with increasing cloud mass M W m- In low- mass halos 
the clouds are pressure-confined, and Mwm is proportional to the cloud volume, so that 

1/3 

r w/H ~ M^ M . The cloud radius is less sensitive to M WM in high-mass halos as the clouds 
become gravitationally confined to approximately fixed volumes. For a given halo mass and 
cloud mass, the cloud radius rw/H decreases as Phim is increased and the clouds are pressure 
confined to smaller volumes. 

The vertically running curve in panels b, c, and f of Figures 6 through 11 identify the 
locus of marginally bound clouds as defined by equation (34). Clouds to the left of this curve 
are unbound and gas particles with thermal energies gas can escape the halo potential wells. 
Clouds to the right of this curve satisfy equation (34) and are fully bound. The position of 
the "bounding" curve shifts towards lower virial masses as Phim increases and the WM gas 
is pressure-confined to the inner parts of the potential wells. For example, 10 4 K clouds in 
median Burkert halos, with Phim = 10 cur 3 K, are bound if M vir £ 1.5 x 1O 8 M (or v s £ 15 
km s _1 ). 

We define twnm as the radius at which the fractional ionization x e = 0.5. This radius 
effectively marks the location of the ionization fronts (when present) and represents the 
transition point from the WIM to the WNM. In critical clouds, r WNM = 0. Critical clouds 
are indicated by the dashed curves in the panels 6c through 11c. For example, for Phim = 10 
cm~ 3 K, and a median 3 x 1O 8 M Burkert halo, a critical cloud forms for Mwm = 2 x 1O 6 M . 
Clouds to the left and below the dashed curves are fully ionized, while clouds to the right 
and above are predominantly neutral. For fixed Mwm, the clouds are ionized when M vir 
is small since the clouds are then extended low-density objects. As M vir is increased the 
clouds are gravitationally confined to smaller volumes, and become denser. Recombination 
becomes more effective, and a neutral core forms. For high-mass halos, neutral cores exist 
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when M WM is large. However, as the WM gas mass is reduced the cloud mass finally becomes 
too small to maintain an optically thick WIM envelope and the cloud becomes fully ionized. 
As Phim is increased the clouds are compressed, and cloud neutralization occurs at smaller 
virial masses. 

We note that gas close to r w / H will be undetectable since by definition the gas columns 
fall to zero at the cloud boundaries. The observed cloud sizes are better represented by the 
gas scale height r gas , or the neutral boundary %nm- Depending on the bounding pressure, 
r gas and twnm can be significantly smaller than rw/H- 111 gravitationally confined clouds 
rwNM is generally a few times larger than the gas scale height, which may be estimated 
using the analytic small- x expressions for r gas (see Appendix A). E.g., for median Burkert 
halos, and assuming c 9i wnm = 8.2 km s" 1 , the analytic gas scale-height ranges from 0.3 to 
0.5 kpc, for M vir between 10 8 and 10 10 M . For -Aa NFW halos, r gas ranges from 0.8 to 0.3 
kpc. Our numerical results are consistent with these analytic estimates. In median Burkert 
halos, and — Aa NFW halos, r WNM ranges from 0.5 to 2 kpc. 

The dependence of the peak HI column density and total HI cloud mass on M vir and 
M W m can now be readily understood. In pressure- confined clouds, the gas densities are 
approximately constant, and the cloud volumes are proportional to the WM gas mass. Thus, 

1 /3 

AHi, P eak ~ M^ M . The pressure-confined clouds tend to be ionized when Phim is low, and 
the HI then exists as a small neutral fraction within the WIM dominated gas. The peak HI 
column is therefore small (< 3 x 10 18 cm" 2 ) and decreases rapidly as M vir decreases, since 
the gas then becomes more highly ionized. In ionized clouds, the total HI mass is a small 
fraction of the total gas mass. For example, for Phim = 1 cm" 3 K, and a M vir = 4 x 1O 8 M 
median Burkert halo with M WM = 8 x 1O 5 M , the neutral HI mass M m = 1O 3 M , or 
~ 0.1% of the total WM mass. Neutral clouds in high- mass halos are in the small-x limit 
with approximately constant volumes. In such clouds the neutral gas dominates the total 
gas mass, and A r H i, P eak is proportional to M WM . 

A sharp increase in A^Hi, pe ak is associated with the transition from ionized to neutral 
clouds as the halo mass increases. For critical clouds, the peak HI column density equals 
3 x 10 18 cm" 2 for both Burkert and NFW halos, independent of M vir or M WM . The fact 
that the peak HI column density in critical clouds equals a nearly constant value of 3 x 10 18 
cm" 2 , independent of M vir , is due to the weak dependence of the critical central densities, 
^•H.ocrit; on the halo mass or scale velocity. This is also indicated analytically (see Table 6 
in Appendix A). For example, for median Burkert halos, a "universal" central gas density 
^H,ocrit ~ 5.5 x 10" 3 cm" 3 is associated with critical clouds, assuming a WIM sound speed 
c 9 ,wim = 11-6 km s" 1 , and J| = 1 (see Table 6). This then implies a fixed ionization 
rate Ccrit = ^H,ocritttB = 1.4 x 10~ 15 s" 1 at the center of critical clouds [where by definition 
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n(e) = n(H)]. Since for our assumed field the unattenuated ionization rate Co = 4 x fCT 14 
s -1 , a well specified opacity is required to reduce the ionization rate to ( crit . The neutral 
fraction then approaches unity for clouds with peak columns exceeding 10 19 cm~ 2 . 

The condition that M^m/Mdm < 1, where Mom is the dark-matter mass within rw/M, 
must be satisfied for self consistent dark-matter dominated solutions. This condition is met 
for the entire parameter space displayed in Figures 6 through 11. The ratio M WM /M DM 
is smallest for low mass clouds in high mass halos. In such objects the gas is restricted 
to the small-x regions where the DM densities are large. A small gas mass then implies 
nnfi/nds «C 1. For fixed M v ; r , the mass ratio increases with Mwm- For a fixed gas mass, 
Mwm/Mqm increases as M v ; r becomes small, since the gas clouds then extend to large radii 
(x ^> 1) where the DM density becomes small. Thus, M WM /M DM is largest in high-mass 
clouds in low-mass halos. The bound clouds are also dark-matter dominated for the range 
of gas masses we consider. Thus, all regions to the right of the bounding curves (indicating 
the marginally bound clouds) in Figures 6 through 11, represent self-consistent, dark-matter 
dominated, hydrostatic solutions for the WM cloud structure. 

In pressure-confined clouds, the central pressure remains small and is insensitive to 
Mwm because Po/Phim is not large in such objects. For example, for low-mass halos, with 
v s ^ 15 km s -1 , Po/k £ 10 cm -3 K, when Phim = 1 cm~ 3 K. In gravitationally confined 
clouds, where the cloud volumes are approximately constant, Po is proportional to Mwm 
and can become very large. E.g., Po/k 10 4 cm -3 K in high-mass halos (v s £ 30 km s _1 ) 
containing Mwm ~5x 1O 6 M . In critical clouds, Po/k ~ 100 cm -3 K, consistent with our 
analytic estimate of nH,ocrit ~5x 10~ 3 cm~ 3 , and equation (26) for the gas pressure. 

We also show (in Figs. 6f-llf displaying the Po/k contours) the regions of parameter 
space where the gas becomes multi-phased at the cloud centers. These regions are indicated 
by the narrow strip defined by the two dashed curves. On the left curve, Po = P m m(0), and 
multi-phased gas is just possible at r = 0. On the right curve, P = P ma x(0), and such clouds 
contain the largest possible multi-phased cores. Since both P min and P max were computed 
assuming heating by the external field, our results reveal the values of M vir and M W m f° r 
which multi-phased behavior is possible in clouds heated by the external field only. It is 
evident that multi-phased cores occur for values of M vir and M W m that are only slightly 
larger than those required for the transition from ionized to neutral clouds. Multi-Phased 
behavior occurs for clouds with peak HI column densities between 2 x 10 19 and 1 x 10 20 
cm -2 . 

To the right of, or above, the multi-phased strip, we have P > P ma x(0). The warm gas 
in the central regions of such clouds must fully convert to CNM, unless additional internal 
heating mechanisms are invoked to alter the local values of P m i n and P max . Thus, for a 
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given halo mass, the Pq = P ma x(0) contour delineates the maximal mass of warm gas that 
can exist as thermally supported WNM heated by the external field. Larger WM masses 
can exist only if internal heat sources drive the gas to the WNM phase. For example, for 
-PhimA = 10 cm -3 K, and a median Burkert halo with M vir = 2 x 10 9 M Q the maximal WM 
mass is 2 x 10 6 M Q (see Fig. 7f). If the gas mass is increased above this value, the heating 
rates at the cloud center are decreased due to the increased cloud opacity, and the gas in 
the core must convert to CNM. This would then drive the total mass of warm gas back to 
the maximum value of 2 x 1O 6 M (leading also to a reduction in the cloud opacity). More 
generally, if the gas becomes multi-phased, with CNM clouds that are not necessarily at the 
center, the total gas mass (including the CNM) can increase beyond the maximum WNM 
mass set by the condition Pq = P max (0). The maximal WM masses increase without limit as 
M vir becomes small, but these clouds are unbound. For example, for bound clouds in median 
Burkert halos, with Pnm/k = 10 cm -3 K, the largest possible WM mass is 1 x 1O 7 M (see 
Fig. 7f). 

The region of the M vir vs. Mwm plane within the "multi-phased strip" , and to the right 
of the "bounding line", represents the region of parameter space within which acceptable 
models for the CHVCs may be found. In this region the WM clouds are bound to gravita- 
tionally dominant dark-matter halos, and contain neutral multi-phased cores. As we discuss 
further below, a key result is our finding that the location of the multi-phased strip coincides 
with the range of peak HI column densities observed in the CHVCs. 



5. Modelling the Dwarf Galaxies Leo A and Sag DIG 

Our main goal is to find halo-cloud models that provide acceptable fits to the HI proper- 
ties of the CHVCs. However, to validate our approach, we first model the pressure-supported 
HI structures observed in the dwarf irregular (dlrr) galaxies Leo A and Sag DIG, both of 
which are well studied members of the Local Group. Modelling the dwarfs is an important 
auxiliary step because their distances are well determined, so that the true physical sizes 
(as opposed to just the angular sizes) of the WNM clouds are known. A specific family of 
halos (e.g. median, —la, —la etc.) can therefore be selected using the observed gas scale 
heights. We list the basic optical properties (Mateo 1998) and 21 cm HI properties (Young 
& Lo 1996; 1997) of Leo A and Sag DIG in Table 1. 

For Leo A we adopt the recently revised distance of 690 ± 60 kpc determined by Tolstoy 
et al. (1998). Young & Lo (1996) assumed a much larger distance of 2.2 Mpc to Leo A based 
on earlier and less certain observations by Hoessel et al. (1994). For Sag DIG we adopt 
Mateo's (1998) recommended distance of 1061 ± 160 kpc (also assumed by by Young & Lo 
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1997). The total V-band luminosities and half-light core radii are 3.0 x 10 7 L and 0.185 kpc 
for Leo A, and 6.9 x 10 7 L and 0.125 kpc for Sag DIG (Mateo 1998). Both galaxies appear 
to have undergone several episodes of star-formation (Tolstoy et al. 1998; Cook 1988). In 
both Leo A and Sag DIG the measured oxygen abundances are O/Hfa 2 x 10~ 5 , implying 
low metallicities Z ps 0.07Z Q (Skillman, Kennicutt & Hodge 1989; Skillman, Terlevich & 
Melnick 1989). 

Young & Lo (1996) carried out high resolution (15") interferometric VLA 21 cm ob- 
servations of Leo A. They found that the optical galaxy is surrounded by a much larger 
2.8 x 1.7 kpc envelope of HI gas. The total VLA 21 cm flux is 65 ± 3 Jy km s -1 , just about 
equal to the single-dish (43-meter Green Bank) flux of 69 ± 3 Jy km s -1 , implying that the 
interferometer detected all of the HI emission. The total HI mass is 7.9 x 10 6 M Q (for d = 690 
kpc). The peak HI column density is 2.7 x 10 21 cm -2 , and emission was detected down to a 
minimum HI column of iVHi,min ~ 2 x 10 19 cm -2 , at a (mean) radius of r(V H i,min) = 1-2 kpc 
from the galaxy center. A key finding of the Young & Lo study is that the global line profile 
can be separated into broad (9.6 ± 0.2 km s _1 ) and narrow (5.2 ± 0.3 km s" 1 ) components. 
The broad (WNM) component contains 80% of the flux and is produced throughout the HI 
image, whereas the narrow (CNM) component is localized near isolated star-forming regions 
within the optical galaxy. In view of these observations we set Af H i, P eak = 2.2 x 10 21 cm -2 
for the WNM distribution in Leo A. 

In their VLA study of Sag DIG, Young & Lo (1997) discovered a similar 3.3 x 3.0 kpc 
HI envelope surrounding the much smaller optical image. The VLA flux is 30.3 ± 1.5 Jy km 
s _1 , close to the single dish flux of 32.6 ± 1.2 Jy km s -1 . The total implied HI mass (for 
d—1.1 Mpc) is 9.3 x 1O 6 M . The HI structure in Sag DIG is somewhat more complex than 
in Leo A, and appears as an asymmetric ring with several clumps of high column density 
gas near the optical galaxy. The peak column density is 1.2 x 10 21 cm~ 2 , and emission was 
detected down to A^min = 5 x 10 18 cm -2 , at r H i,min = 1-6 kpc from the galaxy center. The 
Sag DIG HI line profile can also be decomposed into extended broad WNM (9-10 km s _1 ) 
and localized narrow CNM (4-6 km s _1 ) components. Young & Lo find that 83% of the total 
flux is contained in the WNM component. Thus, A H i iPca k = 1-0 x 10 21 cm" 2 for the WNM 
in Sag DIG 

Given the interferometric HI maps (Fig. 1 in Young & Lo 1996, and Fig. 3 in Young & 
Lo 1997) it is possible to determine the (projected) 1/e gas scale heights for Leo A and Sag 
DIG. In Leo A, averaging over the long and short axes, the projected angular scale height 
is 9 gas ~ 160 arcsec, corresponding to a projected gas scale height of 0.5 kpc for a distance 
of 690 kpc. Determining a gas scale height for Sag DIG is more complicated due to the 
inhomogeneities in the HI emission. However, based on the Young & Lo (1997) map we 



-32 - 



estimate that 6 gas ps 70 arcsec, implying a projected scale height ~ 0.4 kpc. It is significant 
that we find comparable scale heights for Leo A and Sag DIG. 

Can a DM halo model be found, that reproduces (a) the total extents of the HI clouds, 
(b) the observed 1/e gas scale heights, and (c) the peak WNM HI column densities, assuming 
photoionization by the metagalactic field? What is the maximum possible Local Group HIM 
pressure that is consistent with the dwarf galaxy observations? 

For the dwarf galaxies we adopt a representative projected gas scale height of 0.5 kpc, 
and a WNM velocity dispersion of 9 km s -1 . The observed velocity dispersion is consistent 
with ~ 10 4 K gas. Crucially, the observed scale height effectively determines the relative 
concentration of the DM halo, since as is indicated by our analytic expressions for r gas 
(see Table 6, and Figs. 16 and 17, in Appendix A) the intrinsic scale height r gas is weakly 
dependent on the halo scale velocity. Thus, the observed scale height may be used to select 
an appropriate value of a for the dwarf galaxy halos. In particular, it is clear that if the 
dwarfs are embedded in NFW halos, these halos must be significantly underconcentrated, 
with a ?s — 4 relative to the (theoretical) ACDM x vir (M vir ) relation for median halos as 
specified by equation (15). In contrast, median Burkert halos will produce scale heights very 
close to those observed in the dwarfs. It is also evident from our M vir vs. M WM computations 
(Figs. 6-11) for median Burkert halos and —4a NFW halos, that twnm ~ 1 to 2 kpc for HI 
gas column densities ~ 10 21 cm -2 , consistent with the observed cloud sizes and peak HI 
column densities. These results guide the choice of appropriate dwarf DM halo models. 

However, a complication arises in that the optical galaxies do not appear to be dark- 
matter dominated. In Leo A and Sag DIG, (M/Ly) ps 3.5 and 1.4 within the optical radii 
(Mateo 1998). These mass-to-light ratios are within the range expected for evolved stellar 
systems (e.g. Bruzual & Chariot 1993; Sternberg 1998), and are, for example, much smaller 
than observed in the dark-matter dominated dwarf-spheroidals Ursa Minor and Draco where 
(M/L v )q ~ 80 (Aaronson 1983; Gallagher 1994). 

Thus, in modelling the hydrostatic HI gas distributions in Leo A and Sag DIG we 
consider the possibility that stars contribute to the gravitational potential. We make the 
simplifying approximation that the total gravitational potential is the sum 

V? = V-'dm + V?* (35) 

of (independent) DM and stellar potentials. For the stellar potential we adopt a Plummer 
model (see Spitzer 1987). In this model the enclosed stellar mass 

3 

M*(y) = M.. + y = M s J Mjf {y) , (36) 
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where y = r/r s *, and where r s * is the scale radius of the stellar density distribution, and M s * 
is the total stellar mass. We assume that r s * is equal to the observed (projected) half-light 
core radius 14 , and we fix the total stellar mass assuming (M/Ly)© = 1. As representative 
values for Leo A and Sag DIG we set r s * = 0.185 kpc and M s , = 3 x 1O 6 M (Mateo 1998). 

We augment the equation of hydrostatic equilibrium to include the stellar mass. Ex- 
pression (21) for the gas pressure P(x) becomes 



where x = r/r s is the radius in units of the DM halo scale radius r s , and M ds is the halo 
scale mass. 

In Figure 12 we show a specific gravitationally confined model for the dwarfs, in which 
we combine a stellar system (with r s * = 0.185 kpc, and M SJf = 3 x 10 6 M Q ) with a median 
Burkert halo with a scale velocity v s = 30.0 km s _1 corresponding to M vir = 2.0 x 10 9 M Q . 
We assume a low bounding pressure and set Phim = 1 cm -3 K. A good fit is obtained for a 
total gas mass Mwm = 2.2 x 10 7 M Q . Additional model parameters are listed in Table 2. 

In this model the combined stellar and DM mass dominates the enclosed mass at all radii. 
The stellar component (i.e. the optical galaxy) dominates out to the stellar core radius of 0.2 
kpc. At larger radii the DM halo dominates the mass. The WM cloud consists of a WNM 
core with twnm = 1.4 kpc, surrounded by a WIM envelope extending to %/h = 4.4 kpc. 
The WM cloud satisfies equation (34) and is bound. The total HI mass equals 1.3 x 1O 7 M , 
or 61% of the total hydogen mass. Virtually all of the HI mass is contained within r WNM . 
The computed peak HI column density is 1.4 x 10 21 cm" 2 , consistent with the observed peak 
WNM column densities. The computed 1/e scale height for the projected HI column density 
distribution is r Nm = 0.6 kpc, in good agreement with the observed scale heights of 0.5 kpc. 

The HI column falls sharply near r WNM = 1.4 kpc. The HI column decreases from 
2 x 10 19 to 1 x 10 18 cm" 2 , between 1.3 and 1.5 kpc. The sharp decline in the HI column 
densities that occurs between 10 19 and 10 18 cm" 2 is a general feature of the computed gas 
distributions, for both Burkert and NFW potentials, and is due to the photoionization cut- 
off of the HI gas density distribution. The cut-off column of ~ 10 19 cm" 2 is similar to those 
predicted for flattened, rotationally supported, disk systems (e.g. Corbelli & Salpeter 1993; 
Maloney 1993). The radius at which this cut-off occurs is insensitive to the intensity of the 



14 For many stellar density distributions, including the Plummer model, the projected half-mass radius 
r hp = 0.77tVj, where r h is the half-mass radius. For the Plummer distribution r h = 1.3r s *. Thus, 
(Spitzer 1987). 




(37) 
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ionizing radiation field because the ionization front is located at a radius where the total gas 
density is decreasing exponentially. 

For Leo A we are essentially able to reproduce r H i m in = 1-2 kpc for A^ H imin = 2 x 10 19 
cm" 2 as found by Young & Lo in Leo A. In our interpretation this observational cut-off radius 
is very close to twnm, and corresponds to the location of the ionization front. However, if in 
Sag DIG the neutral boundary also occurs near 1.24 kpc, we are unable to account for the 
column of 5 x 10 18 cm" 2 reported by Young & Lo at the significantly larger radius of 1.6 
kpc. In our model we predict a column of only 6 x 10 17 cm" 2 at this radius. As we discuss in 
§6, we encounter a similar difficulty in fitting the extended low column density distributions 
in the CHVCs. 

In Figure 12d we plot the Ha recombination line surface brightness function of 

cloud radius. The Ha emission is produced in the outer WIM envelope, and the computed 
surface brightness includes the equal contributions from both sides of cloud. At the cloud cen- 
ter Iua = 4.3 mR. The surface brightness rises to 8 mR near r WNM due to limb-brightening, 
and then declines rapidly as the projected columns of ionized gas become small. 

The HI mass in our model clearly exceeds the maximum WNM mass, of about5xlO 5 M , 
(see Fig. 6f) that can be maintained by the metagalactic field. This fact is indicated by 
the P m i n and P max curves (shown in Fig. 12e) that were computed assuming heating and 
ionization by the external field. It is evident that P > P max within 0.9 kpc. The central 
pressure Po/k = 6.0 x 10 3 cm" 3 K (corresponding to a gas density n = 0.5 cm" 3 ) is much 
greater than P ma , x (0)/k = 71.5 cm" 3 K. Thus, our model for the WNM in the dwarf galaxies 
implicitly assumes that additional heat sources (presumably associated with the optical 
components) are available to produce the observed mass of warm HI in these objects. We 
note that the computed P m i n and P max are coupled to the ionzation rate which declines with 
increasing cloud opacity. Figure 12f shows that the total hydrogen ionization rate decreases 
from the unattenuated optically thin value of 4.0 x 10" 14 s" 1 at the cloud edge, to 1.0 x 10" 18 
s -1 in the optically thick core. At r = the total ionization rate is 10 times larger than the 
primary rate, due to the dominating effect of secondary ionizations (which are negligible in 
the ionized envelopes). 

The model displayed in Figure 18 is representative. As is indicated by the parameter- 
space study presented in §4, acceptable models for the dwarf galaxy halos can be found 
for other values of M vir , M WM and Prim, and for (underconcentrated) NFW halos rather 
than Burkert halos. It is evident from Figures 6e and 9e that for a low bounding pressure 
{Pmm/k = 1 cm -3 K) acceptable solutions require v s ^ 25 km s" 1 . In such halos, as in our 
particular model, the WNM and WIM components are gravitationally confined. If larger 
bounding pressures are assumed (e.g. Pmu/k = 100 cm" 3 ) lower values of v s are allowed. 
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We note that for our model, the effect of increasing Phim is simply to truncate the outer 
part of the WM cloud at the radius where P gas equals the increased value of Phim- The inner 
part remains largely unaffected by the truncation. For example, if the bounding pressure is 
increased from 1 to f 00 cm" 3 s" 1 , then the WM cloud boundary moves inwards to rw/H ~ 1-3 
kpc. Most of the tenuous WIM mass is removed, but the inner neutral core is preserved. 
However, if Phim is increased above ~ 100 cm -3 K, the neutral core is also truncated, and 
its size shrinks below the observed values. We conclude that the dwarf galaxy HI structures 
are consistent with an upper-limit of Phim A ~ 100 cm" 3 K for the intergalactic medium in 
the Local Group. 

Our model fit demonstrates that the observed extended neutral HI clouds in Leo A 
and Sag DIG may be readily interpreted as dark-matter dominated, pressure-supported gas 
clouds (although we may have difficulty in accounting for extended HI at a level ^5x 10 18 
cm" 2 ). In our model, the central gas density nu,o = 0.46 cm" 3 is much larger than the central 
gas density of nH,ocrit = 5 x 10" 3 cm" 3 in critical clouds, and is sufficiently large to keep the 
gas cloud neutral against the metagalactic radiation field out to ^ 1 kpc. The characterstic 
(or central) DM density in our model is rids = 1-24 cm" 3 (or 3.1 x 1O" 2 M pc" 3 ), in excellent 
agreement with the typical central DM densities of 2.0 x 10" 2 M Q pc" 3 found in HI rotation 
curve studies of rotationally supported dwarfs and low surface brightness galaxies (de Blok & 
McGaugh 1997; Firmani et al. 2000; de Blok, McGaugh & Rubin 2001). Finally, we note that 
while both NFW and Burkert halos may be found that fit the dwarf galaxy observations, the 
Burkert solutions require median halos with "typical" dark-matter concentration parameters 
(in our model x vil = 25.6). However, NFW solutions are problematic in that they must be 
extremely underconcentrated, by ~ — 4<r relative to the ACDM x v i r (M v i r ) relation, to match 
the observed scale heights in the dwarf galaxies. 



6. Modelling the CHVCs 

We have demonstrated that mini-halo models can be constructed for the HI distributions 
in the dwarf galaxies Leo A and Sag DIG. Can similar models be found for the CHVCs? 
Here, we focus on the HI properties of the CHVCs as most recently probed by the high 
resolution single dish observations and interferometric observations presentd by Braun & 
Burton (2000) and Burton et al. (2001) (see also Brims et al. 2001; Putman et al. 2002). 

In their Arecibo observations, Burton et al. (2001) were able to resolve the WNM HI 
gas distributions in a sample of ten CHVCs. They constructed 1° x 1° sized maps of each 
object at a spatial resolution of 3.1 x 3.7 arcmin, and determined the locations of the peak 
HI flux concentrations. They then carried out deep integration "cross-cut" observations to 
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study the detailed gas distributions around the central peaks. Burton et al. found that 
the peak WNM HI column densities range from 2 x 10 19 to 2 x 10 20 cm -2 , and that the 
projected 1/e scale-lengths of the inner parts of the clouds range from 412 to 908 arcsec, 
with a typical value of 690 arcsec. It is significant that the peak HI columns in the CHVCs 
are much smaller than those observed in Leo A and Sag DIG. In their cross-cut observations, 
Burton et al. detected gas down to column densities of ~ 2 x 10 17 cm -2 (corresponding to 
the la sensitivity limit). In some cases low column density "wings" extend out to ~ 1° 
from the cloud centers. In some clouds the HI gas distributions are quite asymmetric, with 
emission at a level of 10 18 cm" 2 persisting out to several gas scale heights on one side of 
the emission peak, with a more rapid decline on the opposite side. The asymmetries could 
be due to the presence of head-tail structures as observed by Bruns et al (2001) in CHVC 
125+41-207, or to an imprecise positioning of the emission centroids. Burton et al. argued 
that intrinsically exponential gas density distributions (but observed in projection) provide a 
reasonable description of the data, although significant departures from this simple density 
law are clearly present. Burton et al. also found that the velocity gradients across the 
clouds are relatively small (~ 10 km s" 1 degree -1 ) implying that rotation is not dynamically 
important in supporting the WNM clouds. 

CNM cores were detected by Braun & Burton (2000) and de Heij et al. (2002c) in 
11 CHVCs (including two that were subsequently observed at Arecibo) in high resolution 
Westerbork interferometric observations. These observations revealed the presence of high 
column density (up to ~ 10 21 cm -2 ) cores displaying narrow (2 — 15 km s -1 FWHM) HI 
line profiles. The narrow linewidth emission clearly traces a CNM component, and typically 
contributes ~ 30% of the total CHVC HI luminosity. Rotation may be important in the 
cores, where velocity gradients are ~ 1 km s -1 arcmin -1 . The radii of the CNM cores range 
from 5 to 17 arcmin, with an average value of 10 arcmin. It is significant that the CNM core 
sizes are comparable to the WNM gas scale heights. 

Given these observational results we adopt the following properties as "typical" for an 
indivdual CHVC, and which should be reproduced in any successful model for the CHVCs; 
(1) a peak WNM HI column of 5 x 10 19 cm -2 , (2) a projected (1/e) WNM gas scale height 
of 690 arcsec for an assumed distance d, (3) HI columns of ~ 1 x 10 18 cm -2 extending to ~ 3 
scale heights, and (4) a multi-phased core region with projected radius equal to or slightly 
smaller than the gas scale height. Properties (l)-(3) are similar to those that fix the dwarf 
halos, with the critical difference that the distances to the CHVCs are not known a priori. 
Property (4) provides an additional constraint, since we wish to construct models in which 
the multi-phased cores arise in gas heated exclusively by the metagalactic field. We will find 
that minihalo models can be constructed that readily reproduce properties (1), (2) and (4), 
while we have difficulty in accounting for property (3). 
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The total number of CHVCs is also an important constraint. We recall that de Heij 
et al. (2002b) identified 67 CHVCs within the Leiden-Dwingeloo survey, and that Putman 
et al. (2002) found an additional 179 CHVCs in their Parkes HI survey, giving a total of 
246 objects. If the CHVCs represent "missing satellites" associated with either the Local 
Group (Blitz et al. 1999), or with the Milky Way (Klypin et al. 1999), then the total 
number of objects should be related to the typical mass of the CHVC mini-halos. The 
numerical simulations of halo evolution predict that the number of embedded dark-matter 
sub-halos contained within any given parent halo increases rapidly with decreasing sub-halo 
mass (Klypin et al. 1999; Moore et al. 1999). As a guide to this predicted behavior we 
adopt the results of Klypin et al. for their ACDM simulation of the Local Group. In their 
computations two massive halos are produced, representing the Milky Way and M31. Each 
massive object contains a large number of embedded sub-halos. Klypin et al. find that the 
number density of sub-halos is proportional to v~ 2 ' 75 , where v s is the sub-halo scale velocity 
15 . They also find that for a given v s , the volume density of objects decreases as ~ 1 / c? 2 , 
where d is the distance from the center of the parent halo. We reexpress the results of Klypin 
et al. (their eq. [3] and [4]) as 

A/-(>^,<rf) = 1.06xl0 3 f M I iT ' p )( / s Y^f-L-) , (38) 
v ' ; V10 12 M / V10 km s" 1 / VIMpc/ ' v ; 

where J\f(> v s ,< d) is the total number of sub-halos with scale velocities greater than v s , 
contained within a distance d from the center of a parent halo with virial mass M viri p. 

With the above observational and theoretical facts in mind, we now consider two basic 
"scenarios" for the CHVCs. We first consider the CHVCs as "circumgalactic" objects asso- 
ciated with the Milky Way halo. In this scenario the characteristic distance to the CHVCs is 
~ 150 kpc. We then consider the possibility that the CHVCs are more distant (d ^ 750 kpc) 
"extragalactic" Local Group objects as proposed by Blitz et al. (1999). As we will show, 
minihalo models for the CHVCs can be constructed for both scenarios, but our analysis 
strongly favors the circumgalactic model. 



6.1. CHVCs as Circumgalactic Minihalos 

Circumgalactic distances for the CHVCs arise naturally if the CHVC halos are drawn 
from the same family of median Burkert halos (or —Acr NFW halos) that successfully re- 
produce the dwarf galaxy observations. If the CHVC halos are not significantly over- or 



15 We note that Klypin ct al. present their results in terms of the maximum circular velocities, w max , of 
the subhalos. We recall that for NFW and Burkert profiles v niax = 0.8v s . 
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under-concentrated with respect to the dwarf halos, the gas scale heights in the CHVCs 
and in the dwarfs must be comparable, since the scale heights are insensitive to the halo 
mass (see the expressions for r gas in Table 6, and Figures 16 and 17, in Appendix A). If 
the scale-heights are in fact comparable, then the CHVCs must be significantly closer to the 
Milky Way than Leo A or Sag DIG. This is because the 1/e angular sizes of the dwarfs (160 
and 70 arcsec) as observed by Young & Lo are significantly smaller than the angular scale 
heights (690 arcsec) of the CHVCs as observed by Burton et al. Given the known distances 
to the dwarfs, the larger CHVC angular scale heights place the CHVCs at distances of 100 
to 200 kpc. At d — 150 kpc, 690 arcsec corresponds to a linear scale of 0.5 kpc. 

The M W m vs. M vir computations displayed in Figures 6 through 11 show that in clouds 
heated by the metagalactic field, multi-phased cores exist for peak HI column between 2x 10 19 
to 1 x 10 20 cm" 2 . The fact that multi-phased cores are predicted for the range of peak WNM 
column densities observed in the CHVCs strongly suggests that the thermal states of the 
CHVCs are controlled by the external field. For a distance of 150 kpc the external field may 
already be dominated by the metagalactic background with a possible (though uncertain) 
contribution from Galactic stellar Lyc photons (see Appendix B). For those models with 
multi-phased cores (i.e., those within the "multi-phased strip" in Figs. 6f through llf) the 
boundaries of the WNM clouds are predicted to occur at r WNM m 0.6 kpc. It is evident 
that many such halo "solutions" for the CHVC behavior can be found for a wide range of 
halo masses and bounding pressures. For example, for Pmm/k = 10 cm -3 K, bound cloud 
models can be constructed for v s > 17 km s" 1 (see Fig. 7) in Burkert halos. Clearly, however, 
gravitationally confined solutions in very massive halos (e.g., setting M vir = 10 10 M Q , with 
M WM = 1O 6 M ) are not relevant, since ~ 200 such objects would exceed the total mass of 
the Milky Way. 

In fact, the large number of CHVCs suggests that they must be pressure-confined in 
low-mass halos. For example, if half of the CHVCs lie within the characteristic distance 
d, then it follows from equation (38) that the CHVC halos must be low-mass objects. In 
particular, for scale velocities v s ~ 12 km s _1 , J\f(> v s ) = 100 within d = 150 kpc, where 
we have assumed M v - 1TtV = 1O 12 M for the Milky Way (see below). Such low scale velocities 
require large confining HIM pressures to keep the gas both neutral and bound to the halos. 

As a specific example of such a pressure- confined solution, we consider a median Burkert 
halo with v s = 12 km s _1 . For this choice M vir = 1 x 1O 8 M . We set Pmu/k = 50 cm" 3 
K, which is the minimum possible bounding pressure that allows 10 4 K WM gas to remain 
bound to such halos. We fix M WM = 1.1 x 1O 6 M in order to produce a peak HI column 
within the range of observed values. For this gas mass, the baryon to DM mass ratio for the 
mini-halos is M WM /M vir = 1.1 x 10~ 2 , and is significantly smaller than the cosmic baryon 
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fraction VL h /Vt m = 0.09. The resulting model structure is displayed in Figure 13, and the 
various model parameters are summarized in Table 3. 

The WM cloud consists of a WNM core with r WNM = 0.67 kpc, surrounded by a 
WIM envelope extending to r w / H = 1.3 kpc. The gas at r w /H is marginally bound with 
W/T=1.0 at the cloud edge (a condition fixed by our choice for Phim)- The central gas 
density nn,o = 2.1 x 10~ 2 cm -3 , and is sufficiently large for the formation of a neutral core. 
The computed HI mass is 2.8 x 1O 5 M , or about 20% of the total WM mass. 

The computed peak HI column is 5.1 x 10 19 cm -2 . The HI column falls to 1/e of the 
central value at r Nm = 0.50 kpc. For an angular 1/e size of 690 arcsec the implied distance 
is d = 150 kpc. The pressure contrast Po/PmM = 4.7, so that Po/k = 235 cm -3 K. The 
central pressure satisfies the condition P m i n (0) < P < P max (0), where P m i n (0)/k = 95 and 
P max (0)/k = 730 cm -3 K. A multi-phased medium is possible out to 0.33 kpc, at which point 
P = Pmm- A multi-phased core of this size is also consistent with the observations (again 
assuming d « 150 kpc). For radii greater than 0.33 kpc, the local gas pressure P falls below 
the local values of P min and the neutral gas must be entirely in the WNM phase. 

We note that for a distance of 150 kpc, clouds in median v s = 12 km s _1 NFW halos 
(as opposed to Burkert halos) would be much smaller than observed. In such systems the 
projected scale height is only 0.2 kpc. In the circumgalactic model, v s = 12 km s _1 NFW 
halos would have to be significantly underconcentrated, with er ps — 2, to give clouds with 
projected scale heights of 0.5 kpc. 

An important prediction of the model displayed in Figure 13 is the presence of a sharp 
drop in the HI column at twnm- As in the case of our dwarf galaxy models, the ionization 
front is located near the radius where the projected column falls to 10 19 cm~ 2 . In our CHVC 
model Ahi = 5 x 10 18 cur 2 at r WNM = 0.7 kpc, and decreases below 1 x 10 18 cm -2 at 
0.8 kpc. We note that in the CHVC model rwNM/^gas = 1-2, whereas in our dwarf galaxy 
model rwNM/^gas = 2.5. This difference is due to the significantly larger HI gas mass, and 
associated peak HI column, in the dwarf compared to the CHVC. 

At this point we encounter a difficulty in accounting for the extended HI gas "wings" 
found by Burton et al. (2001), in at least some of the CHVCs that they observed. To 
illustrate this difficulty, we plot in Figure 13b (as the dashed curve) the projected HI col- 
umn density distribution assuming that the intrinsic HI gas density profile varies as a pure 
exponential out to several scale heights. We recall that Burton et al. suggested that such 
exponentials are representative of the observed gas distributions. We normalize the exponen- 
tial distribution such that the peak HI column density and the projected 1/e scale height are 
equal to those in our CHVC model. In our model, the HI column density decreases to a value 
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of 10 18 cm" 2 at 1.2 (projected) scale heights, whereas for the pure exponential distribution 
the HI column falls to 10 18 cm" 2 at 4 scale heights. Thus, if an exponential density law 
represents the observations down to columns of 10 18 cm -2 , then the photoionization feature 
expected near a column of 10 19 cm -2 would appear to be absent in the data 16 . We caution 
however, that the exponential forms proposed by Burton et al. are only approximate rep- 
resentations. In fact, in some objects (e.g. in the CHVCs 186-31-206 and 230+61+165; see 
their Fig. 11) features that resemble our predicted photoionization cutoffs may be present. 

As is indicated by Figure 13d the predicted Ha surface brightness ranges from 4.3 
mR at the cloud center, to a maximum limb-brightened value of 6.6 mR at a projected 
radius of 0.7 kpc from the cloud center. The predicted Ha surface brightnesses of the 
CHVCs and dwarf galaxies are equal because both cloud systems are optically thick. The 
predicted surface brightness is proportional to our assumed metagalactic field intensity, and 
is consistent with the expectation that in optically thick clouds the emission measure along 
a cloud diameter equals 2nJ*/aB (plus an additional contribution due to the secondary 
ionizations). Direct measurements of the Ha fluxes from the population of CHVCs are largely 
unavailable. However, Weiner et al. (2001) and Tufte et al. (2002) report Ha intensities 
of 20 to 140 mR in five CHVCs. These values are considerably larger than our predicted 
surface brightnesses. We consider it unlikely that the metagalactic Lyc field is more than 
~ 3 times greater than our estimated strength, and additional sources of ionization appear to 
be required. Galactic radiation could contribute if these particular objects are considerably 
closer to the Galaxy than 150 kpc. Collisional ionization at the interfaces between the 
clouds and a dense HIM is another possibility. Further observations are required to determine 
whether the fluxes found by Weiner et al. and Tufte et al. are typical of the entire population 
of CHVCs. 

For v s = 12 km s" 1 halos, neutral bound clouds require large bounding pressures, with 
-PhimA ~ 50 cm" 3 K at distances of ~ 150 kpc from the Galaxy. We now argue that such 
pressures could plausibly be provided by a hot (Thim ~ 2 x 10 6 K) Galactic corona. To 
demonstrate this we construct a simple hydrostatic model for a corona within the Galactic 
potential. We model the Galactic potential as the combination 

V?Gal = V^disk + V?DM (39) 

of a baryonic disk and a dark-matter halo. The coronal HIM pressure at a distance d is then 



16 We remark that purely exponential gas density distributions are expected in NFW halos in the small- a; 
limit (sec Appendix A). However, this requires massive halos. Furthermore, eliminating the photoionization 
cutoff to a column below 10 18 would require reducing our estimate for the metagalactic Lyc flux (for photons 
between 1 and ~ 2 ryd) by a factor of at least ~ 10 3 . 
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given by 

Prim = P mid e [Wal(0) -^ l(d)]/c HiM (40) 
where P mi d is the mid-plane thermal pressure, and chim = 

116T 6 1/2 km s" 1 is the HIM sound 
speed (where T 6 = T/10 6 K). For the disk we assume a cylindrical (R, z) Miyamoto-Nagai 
potential 

^Wirc 

where we set C\ = 8.887 kpc, t> C i rc = 225 km s _1 , a\ = 6.5 kpc, and b\ = 0.26 kpc (Wolfire 
et al. 1995b). In using this expression we fix R to a constant value of 8.5 kpc. At large 
distances from the disk, z/R^> 1, and z ~ d where d is the distance from the Galactic center. 
For the DM halo we adopt a median Burkert halo with M v - a = 1O 12 M . The associated halo 
parameters are x vir = 13.9, r s = 18.6 kpc, n ds = 0.26 amu cm" 3 , and v s = 200 km s~ 2 . The 
virial radius is r vir = 258 kpc, and the circular velocity at the virial radius v vir = 129 km s -1 . 
Our model for the Galactic DM halo is very similar to the model presented by Moore et al. 
(2001). It satisfies the requirements (e.g. Navarro & Steinmetz 2000) that (a) the disk mass, 
Mdisk ~ 10 11 M Q , is consistent with the mass of baryons that could have cooled out of the 
Galactic DM halo, i.e., the requirement that Mdi S k ~ {^b/^m)M vlT (where Qb/^m = 0.09), 
and (b) the dark matter mass within 10 kpc (1.8 x 1O 1O M ) is less than the upper limit of 
~ 3 x 10 10 M Q inferred from Galaxy mass models (Dehnen & Binney 1998). Following Wolfire 
et al. (1995b), we fix the mid-plane pressure such that the coronal HIM produces a total 
X-ray emission measure of 5 x 10" 3 cm" 6 pc (Garmire et al. 1992; Kuntz & Snowden 2000). 
For a coronal HIM temperature of 2 x 10 6 K (corresponding to c 9iH im = 164 km s" 1 ) we find 
that -P m id A = 2-5 x 10 3 cm" 3 K, consistent with the observed mid-plane thermal pressures. 
The total mass of hot coronal gas equals 2.5 x 1O 1O M , and is a significant fraction of the 
baryonic mass of the Galaxy. In Figure 14 we display the resulting coronal pressure as a 
function of the distance d. It is evident that Phim remains large beyond 100 kpc, and ranges 
from 99 to 43 cm" 3 K between 100 and 260 kpc. These pressures are sufficiently large to 
confine and neutralize the gas contained within a population of low-mass minihalos located 
in the circumgalactic environment. 

We emphasize that the coronal HIM pressure maintains the WM gas within the low- 
mass halos in neutral form, and therefore visible as 21 cm sources. If, for example, our 
v s = 12.5 km s" 1 mini-halo containing 1.1 x 10 6 M Q of gas is moved out to larger distances 
where the enveloping HIM pressure is presumably significantly reduced, the cloud expands 
and becomes ionized. Thus, if the CHVCs are infalling objects, a picture emerges in which a 
population of low-mass halos could exist beyond the immediate circumgalactic environment, 
but would be ionized by the metagalactic field and therefore remain invisible as HI sources. 
It is possible that such objects could be detected via absorption line studies of ionized gas 
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tracers (e.g. Sembach et al. 1999). As such objects approach the higher pressure environment 
of the Galaxy the clouds are compressed and neutralized, and become visible HI sources. 

A distance of d — 150 kpc is likely well within the virial radius of Galactic DM halo 
(for our Galaxy model, r vir = 260 kpc). Our model CHVC mini-halo must therefore be 
tidally truncated at a radius that is smaller than the (original) mini-halo virial radius. What 
is the tidal truncation radius at 150 kpc? Assuming the CHVC is on a circular orbit at 
a distance d around the Milky Way, the truncation radius r t is set by the tidal criterion 
3nM\v(^) — ^CHVc( r t) where nuw{d) is the mean density of the Galaxy within d, and 
ncHvc( r t) is the mean DM density of the CHVC halo within r t . For our Galaxy model, 
3%w — 2.4 x 10~ 3 cm -3 at 150 kpc. For our v s = 12 km s" 1 minihalo, the implied tidal 
radius r t = 6.8 kpc, well within the mini-halo virial radius of 12.0 kpc. However, r t is still 
significantly larger than the WM cloud radius r w / H = 1.3 kpc set by the confining HIM 
pressure. The actual truncation radius will depend on the orbital parameters, and could be 
significantly smaller than our estimate if the periastron distances are sufficiently small. If 
the CHVCs come too close to the Galaxy, they will be tidally disrupted. As we now discuss, 
they may also be destroyed by ram-pressure stripping. 

While low-mass halos require a large HIM pressure to keep the gas bound and neutral, 
if the HIM becomes too dense the CHVC gas will be ram-pressure stripped out of the 
mini-halos. Stripping occurs when the ram-pressure exceeds the gravitational force per unit 
area on the gas (Gunn & Gott 1972; Blitz & Robishaw 2000; Mori & Burkert 2000; Quilis 
& Moore 2001). For neutral clouds moving with orbital velocity v through a hot ionized 
medium, stripping will occur when 

( nmM ) >4A( 5™ )( Nm )( - Y 2 (42) 

Vl0- 4 cm- 3 y~ ' \10 7 M o kpc -2 / V10 20 cm-V V100 km s- 1 / 1 } 

where riuiM is the HIM gas density, Edm is the dark-matter surface density, and Ahi is the 
HI column density. For our model CHVC the DM surface density within the projected 1/e 
scale height is Edm = 1-0 x 10 7 M Q kpc~ 2 . Setting A%i equal to the peak column density 
of 5.1 x 10 19 cm -2 , and assuming a typical CHVC velocity of 100 km s -1 , it follows that 
stripping will occur for izhim ~2 x 10 -4 cm~ 3 . This is much larger than a HIM density of 
2.5 x 10 -5 cm~ 3 for a pressure Pmu/k = 50 cm~ 3 K, and Xhim = 2 x 10 6 K. We conclude 
that our pressure-confined v s = 12 km s™ 1 mini-halo model clouds (and the cores of the 
observed clouds) would be stable against rapid stripping by our postulated coronal HIM at 
150 kpc. We note, however, that at HIM densities of 2.5 x 10~ 5 cm -3 , the observed (but not 
predicted) low-column density (^5x 10 18 cm~ 2 ) wings, should be stripped off the clouds. 
Their presence thus poses a challenge to our view of the CHVCs as pressure-confined objects. 
The cloud wings could conceivably represent stripped gas, but then the question remains as 
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to why this low column density gas is not photoionized by the background field (and by the 
additional sources of ionization indicated by the Weiner et al. observations). 

We mention one additional problem. If the bounding pressure is reduced sufficiently, 
not only do the clouds become ionized, they also become unbound. For example, it is clear 
from Figure 6c that for Pmu/k = 1 cm~ 3 K, our CHVC model lies in the unbound portion of 
the Mwm vs M vir diagram, indicating that mass-loss from the outer layers of the cloud would 
occur for such a low HIM pressure. Thus, the question arises as to how the gas was retained 
by the mini-halos prior to their approach into the circumgalactic environment. This problem 
can be alleviated if the graviational potential well is deepened by increasing the mini-halo 
mass. E.g., increasing v s to 17 km s _1 is sufficient to keep our WM mass of 1.1 x 1O 6 M 
bound (although still ionized) for Pum/k = 1 cm -3 K (see Fig. 6c). However, the expected 
number of objects decreases rapidly as v s is increased as indicated by equation (38). 



6.2. CHVCs as Extragalactic Minihalos 

We now consider the possibility that the CHVCs are more distant "extragalactic" ob- 
jects, and represent sub-structure on the scale of the entire Local Group (LG). This is the 
basic picture proposed by Blitz et al., who placed the CHVCs at characteristic distances 
of >, 750 kpc. As we will show, models for the CHVCs can be constructed for such dis- 
tances. However, there are two major difficulties with this picture. First, the mini-halos 
must be extremely underconcentrated. Second, the reduced central gas pressures suppress 
the formation of multi-phased cores. 

For d = 750 kpc, the observed CHVC HI gas scale height of 690 arcsec corresponds 
to a linear scale of 2.5 kpc. Such large values for the gas scale height require extremely 
underconcentrated halos, for both Burkert and NFW potentials (see Table 6 and Figs. 16 and 
17). Such objects must deviate very substantially from the correlation x v i r (M vir ) expected 
for median ACDM halos, with o £ — 4 for Burkert halos, and o ^ — 6 for NFW halos. 
This is a major problem with the LG hypothesis. At distances ^ 750 kpc, the CHVC halos 
would have to be quite different from those found in the CDM collapse simulations (see §2.2). 
From a more observational perspective, the CHVCs would then represent objects that are 
quite distinct from the dwarf galaxies. For example, for r gas = 2.5 kpc the required DM scale 
density in Burkert halos (see Table 6) must be extremely small, with n ds w 5 x 10~ 2 cm~ 3 . 
This is much smaller than the characteristic DM densities of ~ 1 cm -3 that we find for Leo 
A and Sag DIG, or that have been inferred from the dwarf galaxy rotation curve studies. 

As a specific example of a CHVC model for d ps 750 kpc, we again consider a pressure- 
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confined system. We select a — Aa Burkert halo with v s = 12 km s" 1 . The corresponding 
virial mass is 3 x 1O 8 M . We set the confining pressure equal to the minimum possible HIM 
pressure for bound cloud solutions. For this model the minimum pressure is Pmu/k — 20 
cm -3 K. We choose a WM gas mass equal to Mwm — 5.8 x 1O 7 M , to give a total HI column 
within the observed range. The resulting model structure is displayed in Figure 15, and the 
model parameters are summarized in Table 4. We note that for this model the concentration 
parameter x vir = 8.1. This is extremely small for our assumed halo mass. For median halos, 
such concentrations are expected for masses ~ 1O 13 M (see eq. [14]), corresponding to the 
mass scale of an entire galaxy group. 

In our extragalactic model, the peak HI column density equals 1.0 x 10 20 cm -2 . The 
intrinsic WM gas scale height is r gas = 3.0 kpc, and the projected 1/e HI gas scale height 
r Nm = 2.6 kpc. For a WNM core size of 690 arcsec, the implied distance is 777 kpc. The 
pressure contrast Po/-Phim = 4.4, so that Po/k = 88 cur 3 K. The central pressure satisfies 
the condition for multi-phased behavior P m i n (0) < P < P max (0) where P m i n (0) = 67 and 
-fmax(O) = 523 cm" 3 K. We note that due to the larger distance, the central gas pressure 
is significantly lower than in our circumgalactic model. Large peak HI column densities are 
therefore required to provide the opacity necessary to reduce P m i n and P max to values that 
allow multi-phased gas in the core. For —4a Burkert halos we find that multi-phased gas is 
possible for columns between 9 x 10 19 to 3 x 10 20 cur 2 . In our specific extragalactic model a 
small multiphase core is just possible out to 1.2 kpc, consistent with the observed presence 
of CNM cores within the WNM gas scale height. While the required WNM column of ~ 10 20 
cm" 2 is within the range of observed values, most of the CHVCs are not as shielded as this 
(Burton et al. 2001). At large distances, the multi-phased cores in most objects would 
vanish. 

The WNM component extends to r WNM = 3.4 kpc, and is surrounded by a WIM envelope 
that extends to 7.2 kpc. The total HI mass is 1.4 x 10 7 M Q . The central gas density n H ,o = 
7.6 x 10" 3 cm -3 is sufficiently large for the formation of a neutral core (since for r gas m 2.5 
kpc nn.ocrit = 2 x 10" 3 cm" 3 ). As expected, the ionization front is associated with a projected 
HI column close to 10 19 cm" 2 . At 3.0 kpc from the cloud center, Nm = 8 x 10 18 cm -2 and 
decreases below 10 18 cm -2 at 4.0 kpc. Extended low column density wings at a level of 
~ 10 18 cm" 2 are not present. 

If the CHVCs are at large distances, and they are pressure-confined objects, their large 
extents set limits on the intergalactic medium (IGM) pressure within the Local Group. The 
bounding pressure Phim cannot exceed ~ 40 cm" 3 K, since otherwise the clouds would be 
compressed to radii less than 2.6 kpc. Conversely, for v s = 12 km s _1 the bounding pressure 
cannot be smaller than 20 cm -3 K, since otherwise the clouds would become unbound. The 
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lower limit on the pressure can be reduced if more massive halos are considered. If the IGM 
pressure is close to the required ~ 20 cm -3 K this gas must also be very hot. For a uniform 
density HIM contained within a sphere of radius d the total HIM mass is 

M HIM ^,0x 1 0"( r ^)( I ^_)(^|)- 1 M (43, 

where T mM is the HIM temperature. Assuming the HIM mass does not exceed a fraction 
~ 0.25 of the total baryonic mass ~ 3 x 1O 11 M of the LG (essentially the Milky Way plus 
M31) then for d = 750 kpc, and PhimA = 20 cnT 3 K, it follows that T fflM ~ 10 7 K. It is 
difficult to see how such hot gas could remain bound to the Local Group. 

Finally, another important feature of models for d ^ 750 kpc is a much larger WM to 
DM mass ratio. In our extragalactic model, M WM /M vir = 0.22 (more than twice the cosmic 
baryon fraction) as opposed to 0.011 in our circumgalactic model. For v s = 12 km s -1 , M vir 
grows by a factor of 3.2 as a is reduced from to -4. In contrast, M WM ~ nmnNmr^, so 
that Mwm grows by a factor £ 50. The gas-to-DM mass fraction can be reduced if v s is 
increased. However, it is clear that as the CHVC halos are placed at larger distances they 
must become increasingly gas rich. This poses an additional problem for the extragalactic 
hypothesis, as it is difficult for low-mass halos to retain their gas. 

7. Summary 

In this paper we examine the hypothesis that the compact high velocity HI clouds 
(CHVCs) trace sub-structure within the dark-matter halo of the Galaxy, or within the entire 
Local Group system. 

For this purpose we carry out detailed computations of the coupled hydrostatic and ion- 
ization structures of pressure supported clouds that are confined by gravitationally dominant 
dark-matter (DM) "mini-halos" . We focus on low-metallicity systems that are ionized and 
heated by the metagalactic background field. We provide a fit for this field (from the near-IR 
to the hard X-ray regime) that is based on a combination of observations and theoretical 
estimates. We explicitly consider the effects of external bounding pressures on the mini-halo 
cloud structure. We consider bounding pressures provided by a hot ionized Galactic corona 
or an intergalactic medium in the Local Group. 

We consider dark-matter halos with either cuspy (NFW) or constant density (Burkert) 
cores. We adopt a ACDM cosmological model to parameterize the expected correlations 
between the various halo scale parameters, as well as the dispersion in the correlations. We 
focus on low-mass halos with virial masses ~ 10 8 to 10 10 M o , (or halo scale velocities between 
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~ 10 to 50 km s 1 ) containing warm (10 4 K) clouds of neutral plus photoionized gas with 
masses between ~ 10 4 and 1O 8 M . 

We determine how the cloud sizes and hydrogen gas distributions, as well as the gas 
phase-states - ionized (WIM), neutral (WNM), or multi-phased (WNM/CNM) - depend on 
the halo parameters, gas masses, and bounding pressures. We consider both gravitationally 
confined systems, where the gas is effectively restricted to within the halo scale radii (or 
halo cores), as well as pressure-confined clouds in which the gas can extend to large halo 
radii. We find that gravitationally confined systems are those for which v s /c g ^ 1.5, where 
v s is the halo scale velocity and c g is the gas sound speed. We determine the conditions 
required for the transition from ionized to neutral clouds, and for the formation of thermally 
unstable multi-phased (WNM/CNM) cores within neutral clouds. We also compute the 
emission measures and Ha recombination line surface brightnesses that are produced in the 
WIM components of the halo clouds. For optically thick clouds ionized by the metagalactic 
field, we find that the Ha surface brightness ranges from 4 mR at the cloud centers to 7 mR 
at the limb-brightened edges. 

As a step in our computations we construct pressure vs. density phase diagrams for 
low metallicity (Z = 0.1Z Q ) gas that is heated by the metagalactic field for a wide range 
(10 18 — 10 21 cm" 2 ) of assumed HI shielding columns, and we compute the critical phase- 
transition pressures P m i n and P max as functions of the hydrogen ionization rate. We use the 
results of these computations to determine the thermal phase states of the neutral gas in our 
halo cloud models. 

We present mini-halo models for the pressure supported HI structures observed in the 
Local Group dwarf irregular galaxies Leo A and Sag DIG. We include the stellar contribution 
to the gravitational potentials. We identify the HI cloud boundaries observed in Leo A and 
Sag DIG with the ionization fronts, and we derive an upper limit of Phim/ k ^ 100 cm -3 K for 
the ambient pressure of the intergalactic medium in the Local Group. The distances to Leo 
A and Sag DIG are well established, and the observed HI gas scale heights of 0.5 kpc in these 
objects imply characteristic DM densities of 1.2 amu cm" 3 (or O.O3M pc" 3 ) for the DM 
halos. These densities are consistent with those previously found via rotation curve studies 
of rotationally supported dwarfs and low-surface brightness galaxies. Leo A and Sag DIG 
obey the halo correlations that are expected for typical ("median") DM halos in a ACDM 
cosmology, provided the halos contain constant density cores. NFW halos would have to be 
extremely underconcentrated (by — 4<r relative to median halos) given the observed gas scale 
heights. 

For the family of median Burkert (or — 4er NFW) halos that obey the scalings and 
correlations as fixed by the dwarf galaxies, we find that the transition from ionized to neutral 
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clouds occurs for central WNM HI column densities equal to 3 x 10 18 cm -2 , independent of 
the halo virial mass or total gas mass. The gas is fully neutral within projected radii where 
HI column exceeds 10 19 cm -2 . A photoionization cut-off is predicted near 10 19 cm -2 (similar 
to what is found for flattened, rotationally supported, disk systems) independent of the halo 
mass. We find that for clouds heated by the metagalactic field, multi-phased cores occur for 
peak WNM columns between 2 x 10 19 and 2 x 10 20 cm' 2 . The peak WNM columns of ~ 10 21 
cm~ 2 observed in the dwarf galaxies require additional (possibly internal) heat sources. 

We construct explicit "mini-halo" models for the multi-phased (and low-metallicity) 
compact high-velocity HI clouds. We consider the CHVCs as either "circumgalactic" objects 
associated with the Milky Way halo, with implied characteristic distances of 150 kpc, or as 
more distant £ 750 kpc "extragalactic" Local Group objects. If the CHVC halos are drawn 
from the same family of halos that successfully reproduce the dwarf galaxy observations, 
then the CHVCs must be circumgalactic objects. The observed 1/e gas scale heights then 
correspond to sizes of 0.5 kpc (as in the dwarfs). The typical HI mass for an individual 
CHVC is then ~ 3 x 10 5 M Q , and the total HI mass in the entire population of CHVCs is 
~ 6 x 1O 7 M . The predicted multi-phased behavior that occurs for peak WNM columns 
between 2 x 10 19 and 2 x 10 20 cm -2 is consistent with the observed multi-phased behavior and 
range of peak WNM columns in the CHVCs. If the large population of CHVCs represent 
"missing low-mass DM satellites" of the Galaxy, then these HI clouds must be pressure- 
confined to keep the gas neutral within the weak DM potentials. The observed number of 
CHVCs imply typical CHVC mini-halo scale velocities of v s = 12 km s -1 . For such objects 
the confining pressure must exceed ~ 50 cm' 3 K to keep the gas bound. We construct a 
simple model for a Galactic corona, and show that a hot 2 x 10 6 K Galactic corona (in 
pressure equilibrium with the Galactic disk) could provide the required pressure at 150 kpc. 

If the CHVCs are "extragalactic" objects with distances £ 750 kpc, then their asso- 
ciated halos must be very "underconcentrated" , with characteristic DM densities ,$,0.08 
cm -3 , much smaller than expected for their mass, and significantly smaller than observed in 
the dwarf galaxies. Multi-Phased cores are possible, but require shielding columns that are 
generally higher than observed. For d = 750 kpc, the typical CHVC HI mass is ~ 1 x 1O 7 M , 
and the total HI mass in the system of CHVCs is ~ 2 x 1O 9 M . 

Our analysis favors the circumgalactic hypothesis for the location of the CHVCs. In 
this picture, the CHVCs represent pressure-confined clouds that are associated with tidally 
truncated dark-matter sub-halos that have survived the hierarchical formation process of the 
Galactic halo. The CHVC mini-halos are visible as neutral 21 cm sources due the compression 
provided by an ambient Galactic corona. 

Our analysis also appears to favor Burkert halos with constant density dark-matter 
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cores, as opposed to NFW halos with diverging central DM densities. While Burkert halos 
are consistent with the observed scale heights at circumgalactic distances of 150 kpc, NFW 
halos would have to be significantly underconcentrated to yield the observed cloud sizes for 
such distances. 

We compute the maximum masses of WNM gas that can be maintained by the meta- 
galactic field within the DM halos. Our analysis suggests that in the mini-halo scenario 
for the CHVCs, these objects are "failed galaxies" and did not form stars, because most 
of the neutral gas was always maintained as WNM by the metagalactic field. In contrast 
the halos associated the dwarf galaxies in the Local Group (and elsewhere) contained suffi- 
ciently large masses of gas such that most of the neutral component was converted to CNM, 
a precondition for star-formation. 

Our hydrostatic mini-halo cloud models are able to account for many properties of the 
CHVCs, including their observed peak HI columns, core sizes, and multi-phased behavior. 
However, important theoretical and observational difficulties remain. Theoretically, a ques- 
tion remains as to the origin of the gas in the CHVC mini-halos. The gas is unlikely to 
be purely primordial, since such gas is easily lost from low-mass mini-halos as the bound- 
ing pressure becomes low. The presence in some objects of extended low column density 
HI wings, and Ha emission line fluxes in several CHVCs that are significantly larger than 
expected, represent significant observational challenges to the mini-halo models we have 
presented in this paper. Additional high resolution HI mapping observations and sensitive 
Ha line measurements are required to determine the possible contributions of gas stripping 
and perhaps collisional ionization to the CHVC cloud structure. Such observations will help 
establish, or possibly refute, the hypothesis that the CHVCs trace dark-matter substructure. 

We thank Leo Blitz, Stu Bowyer, James Bullock, Butler Burton, Marc Davis, Orly Gnat, 
David Hollenbach, Tsafrir Kollat, and David Spergel for many helpful discussions during the 
course of this work. We thank James Bullock for providing us with his CDM halo code. We 
also thank the referee for comments and suggestions that improved this paper. C.F.M. was 
supported in part by NSF grant AST-0098365. M.G.W. was supported in part by a NASA 
LTSA grant NAG5-9271, and by NSF grant AST 95-29167. 

A. NFW AND BURKERT HALOS AND THE SMALL-X LIMIT 

In this appendix we derive some useful analytic formulae for the cloud properties and 
gas distributions when the clouds are confined to the "small-x" (x < 1) regions of the 
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dark-matter halos. 

We first summarize the basic properties of the NFW (Navarro et al. 1997) and Burkert 
(1995) dark-matter halos. In Table 5 we list the expressions for the (dimensionless) density 
profile, f p , the enclosed mass, /m, and the gravitational potential, f v , (as defined in §2.1) 
for NFW and Burkert halos. 

In NFW halos the central densities diverge, with f p ps l/x for x <C 1, while the enclosed 
mass and potential remain finite, with f M (3/2)x 2 and f v ps (3/2)x. The scale density p ds 
(see §2.1) equals the local DM density at x = 0.46. At x — 1 the DM density pd = Pds/4, 
and the enclosed mass Md = 0.58Md s . The enclosed mass Md = Md s at x = 1.6. For x 1, 
/ p ~ l/x 3 , ju diverges logarithmically, and f v approaches an asymptotic value of 3. The 
circular velocity v = v s f v reaches a maximum value t> max = 0.8t> s at x = 2.2. 

Burkert halos are defined such that they contain constant density dark-matter cores. 
Thus, for x <C 1, f p ~ 1, /m ~ x 3 , and f v ~ x 2 /2. The scale density pd s equals the 
central (core) density Pd(0). At x — 1 the DM density pd = Pds/^, and the enclosed mass 
Md = 0.38Md s . The enclosed mass Md = Md s at x — 1.7. For x ^> 1, f p ~ l/x 3 , diverges 
logarithmically, and f v approaches an asymptotic value of 37r/4. The circular velocity v = 
v s f v reaches a maximum value f max = 0.8v s at x = 3.3. 

We now consider simplified isothermal clouds with fixed ionization state (e.g., fully ion- 
ized or completely neutral) so that the gas is characterized by a constant velocity dispersion 
c g . Writing equation (20) as 



it follows that the hydrostatic gas density profile for a dark-matter dominated cloud is 



where f v is the dimensionless halo potential. For a given halo profile (NFW or Burkert) the 
gas density distribution depends on the single parameter v s /c g . In Table Al we include the 
resulting analytic solutions for the gas distributions (Makino, Sasaki & Suto 1998; Wu 2000, 
Wu & Xue 2000) as given by equation (A2). 

For a; C 1 the expressions for the gas distributions simplify considerably, and simple 
analytic formulae for the gas scale heights, column densities, masses, and cloud emission 
measures may be derived. We list these various formulae in Table 6. 

In the small-x limit, the gas density distribution / gas varies exponentially in NFW halos. 
In Burkert halos the gas distribution is a Gaussian. For both halo types, we define the gas 
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scale height, x gas = r gas /r s , as the radius at which the gas density decreases to a factor 
1/e times the central value. Expressions for r gas as functions of the DM halo parameters 
are listed in Table 6. In NFW halos r gas depends on two (correlated) halo scale paramters, 
e.g. rids and r s , or r s and v s . In Burkert halos, r gas may be expressed as a function of 
only, that is, on the value of the uniform dark-matter density in the i C 1 halo core. 

In Figures 16 and 17 we plot the small- x values of r gas as a function of the halo scale 
velocity, v s , for median, +3a over-concentrated, and — 3cr under-concentrated NFW and 
Burkert halos. In Figures 16 and 17 we display r gas as a function of M vir . In these plots 
we assume a (WNM) sound speed of 8.2 km s" 1 . Analytic expressions for r gas as functions 
of a and v s , or as functions of a and M w - U may be obtained using the scaling relations (18), 
(19), and (17) for n ds , r s and M vir (see §2.2). We list these expressions for r gas in Table 6. 
In both NFW and Burkert halos the scale heights for median (a = 0) halos are small. For 
median NFW halos r gas ranges from 0.1 to 0.04 kpc for v s between 10 to 50 km s -1 , or virial 
masses between 10 8 and 1O 1O M . For median Burkert halos r gas ranges from 0.2 to 0.6 kpc 
for v s from 10 to 50 km s _1 , or M v ; r between 10 8 and 1O 1O M . We note that for Burkert 
halos r gas is a very weak function of M v i r . The scale heights may be increased by moving 
to under-concentrated halos. We illustrate this in Figures 16c, 16d, 17c, and 17d where we 
display level curves for r gas in the range of 0.5 to 2.5 kpc, in a vs. v s , and a vs. M vir diagrams. 
It is apparent that large scale heights require extremely underconcentrated halos. 

In the small-x limit the total gas mass remains finite as £w/h becomes arbitrarily large 
(i.e. as the bounding pressure is reduced to zero). Similarly, the gas column densities at 
all impact parameters b = r p /r gas (where r v is the projected radius) remain finite. Small-rr 
expressions for the cloud masses and columns are listed in Table 6. In NFW halos the column 
densities are proportional to bK 1 (b) where K 1 is the first-order modified Bessel function. The 
product bK 1 (0) = 1 for b = 0, and declines by a factor of 1/e at b = 1.65. Thus, for NFW 
halos the projected 1/e scale height is 1.65 times larger than the intrinsic gas scale height. 
In Burkert halos the column density distribution is a Gaussian, and the projected 1/e scale 
height is equal to the intrinsic scale height. 

Analytic conditions for the transition from ionized to neutral clouds may be derived 
based on the behavior of the cloud emission measure with varying optical depth. The emis- 
sion measure depends on the fraction of the incident flux, ttJ*, of Lyman continuum photons 
that are absorbed by the cloud. In optically thick clouds, the entire Lyc flux is absorbed 
in a thin shell, and EM(x = 0) = 2nJ*/aB- In the opposite limit, EM must vanish in 
optically thin clouds as the cloud mass becomes sufficiently small. Our radiative transfer 
computations show that EM always reaches a maximum value of 
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in critical clouds where the fractional ionization x e = 0.5 at the cloud center. Indeed, 
it is straightforward to show analytically that for spherical clouds (and with the neglect 
of secondary ionizations) the emission measure must satisfy the upper and lower bounds, 
87iJ*/aB > EM > 2nJ*f/aB, where / is the fraction of the incident radiation that is 
absorbed. Our radiative transfer computations are consistent with these limits. We find 
that EM vanishes when the cloud is optically thin and / is small. The emission measure 
then rises to a maximum value of close to An J* /a B in critical clouds, and then declines to 
EM = 2-kJ* /a B in optically thick (/ = 1) clouds. 

The emission measures for fully ionized clouds are listed in Table 6. Crucially, these 
emission measures cannot exceed i?M max . It follows that to remain fully ionized the gas 
scale height must satisfy an upper bound, with 



where / = 1 for NFW halos, / = \/2Jti for Burkert halos, and J| = J*/10 3 photons cm' 2 
s _1 sr _1 (for our assumed field J| = 1.02). Neutral cores must form when r gas exceeds this 
limit. We identify critical clouds as those for which equality holds in equation (A4). Given 
our formulae for r gas equation (A4) implies a critical central gas density, nH,ocrit above which 
neutral cores must form. The resulting expressions for nH,ocrit, as functions of v s or M v ; r , are 
listed in Table 6. We note that in evaluating these critical densities, a WIM sound speed 
should be adopted since critical clouds are WIM dominated. It is evident that for NFW 
halos the critical density is a weakly increasing function of the halo scale velocity or mass, 
whereas in Burkert halos the critical density is a very weakly decreasing function of the scale 
velocity or mass. 

For what values of v s /c g is the small-x approximation valid? It follows from the analytic 
expressions for / gas listed in Table 5 that for v s /c g ^ 2 the gas densities become vanishingly 
small for x ^ 1, while for v s /c g £ 2 the gas densities remain relatively large out to large halo 
radii. Given the analytic expressions for / gas listed in Table 5 we have numerically integrated 
the total gas mass, MwM/47rr 3 nH,o, peak gas column, A 7 wM P eak/2nH i o r s, and emission mea- 
sure EM/2n^ r s (assuming fully ionized clouds) as functions of v s /c g . We then compared 
the results of the numerical integrations to the values as given by the small-rr formulae. We 
find that for v s /c g > 1.5 the small-x expressions for the gas column (and therefore also the 
gas scale height r gas ), and the emission measure, are good to better than a factor of 2, for 
both NFW and Burkert profiles. For v s /c g > 2 the accuracy is better than a few percent. 
The total gas mass diverges rapidly at small v s /c g , and the small-x expressions for the total 
gas mass become accurate only for v s /c g > 3. For such large values of v s /c g the gas mass is 
effectively independent of the bounding pressure as Phim becomes arbitrarily small. 
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B. METAGALACTIC RADIATION FIELD 

In this paper we construct models assuming the (z = 0) metagalactic field is the domi- 
nant source of radiation impinging on the mini-halo clouds. 

The external field, J u (ergs cm -2 s" 1 Hz -1 sr" 1 , hereafter cgs units), is a critical quantity 
controlling the overall cloud structure. For a given total gas mass, the Lyc flux determines the 
fractions contained in the WIM vs. the shielded WNM components. The radiation field also 
determines the thermal phase state of the HI gas (see §3.2), since hard ultraviolet (EUV) and 
X-ray photons can penetrate into the shielded cores and heat the HI gas via photoionization. 
Non-ionizing radiation is also important. Far-ultraviolet (FUV) grain photoelectric emission 
is a heat source, while FUV photons produce the C + ions which dominate the cooling of the 
CNM. Molecular hydrogen (H 2 ) is photodissociated by the FUV field, and the abundance of 
H~, a key intermediary in the gas phase production of H 2 , is limited by optical and infrared 
photodetachment. 

In view of the above processes the metagalactic field must be specified from X-ray 
energies ~ 1 keV, above which clouds with the largest HI columns we consider (~ 10 22 cm~ 2 ) 
become optically thin, down to to the near- infrared H photodetachment threshold at 1.64 
fim. We adopt a representation for the metagalactic field that is based on a combination 
of observational constraints and theoretical estimates. We then estimate the distance from 
the Galactic disk at which the metagalactic field dominates. We first consider the ionizing 
radiation, and then the FUV through IR components. 

The metagalactic Lyc flux is generally believed to be produced by quasars with an 
additional possible contribution from star-forming galaxies. The background field is directly 
observable at X-ray wavelengths (Gendreau et al. 1995; Chen, Fabian & Gendreau 1997; 
Miyaji et al. 1998; Hasinger 2000), but is unobservable near the Lyman limit due to efficient 
absorption by Galactic HI. 

The unobservable flux at z = can be constrained via Ha recombination line measure- 
ments from (starless) intergalactic HI clouds (Vogel et al. 1995; Weiner et al. 2001), or via 
measurements of the HI "truncation column densities" at the edges of galaxy disks (Sunyaev 
1969; Silk 1971; Maloney 1993). Vogel et al. (1995) determined a (2a) upper limit of 20 mR 
for the Ha surface brightness in the intergalactic Giovanelli & Haynes cloud. This measure- 
ment sets an upper limit on Arc J* (eq.[25]). For opaque spherical clouds, An J* = a B {EM) 
where (EM) is the emission measure averaged over the projected surface area. We conclude 
that An J* < 4.5 x 10 4 photons cm~ 2 s _1 in the local IGM. 

Theoretical calculations of the metagalactic ionizing radiation field at low redshifts have 
been presented by Haardt & Madau (1996), Okoshi & Ikeuchi (1996), and Shull et al. (1999) 
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with similar results. At z = Haardt & Madau estimate that J u0 = 2.0 x 1CT 23 cgs at the 
Lyman limit, while Shull et al. find J u0 = 1.3 x 10~ 23 cgs. The Lyc intensity drops sharply 
(as an approximate power-law) to a value of ~ 3 x 10~ 25 cgs at the Hell edge (4 ryd), and 
flattens at higher energies where the IGM becomes less opaque. 

In estimating the background field below 4 ryd to the Lyman limit we rely on on the 
Haardt & Madau calculations. At energies greater than 4 ryd including EUV and X-ray 
wavelengths we rely on observations. Chen et al. (1997) (see also Gendreau et al. 1995) find 
that between 1 and 7 keV the X-ray background (as observed by ROSAT and ASCA) is well 
fit by a simple power-law, AE^, where Ek ev is the photon energy in kev, A = 10.5 ± 0.3 
photons s -1 cm -2 kev -1 sr _1 , and the photon index T = 1.46 ± 0.06. A very similar result 
was found by Miyaji et al. (1998). At soft X-ray energies ( ^ IkeV) it becomes difficult to 
disentangle the background from thermal emission produced in the (distant) Galactic halo 
and (nearby) Local Bubble of hot gas. However, in an analysis of soft X-ray "shadows" 
(i.e. anticorrelations on the sky between soft X-ray flux and absorbing HI column density) 
Snowden et al. (2000) argue that the Chen et al. power-law plausibly extends down to 1/4 
keV, where additional flux (of distant origin) at 1/4 keV is due to a thermal 10 6 ' 4 K Galactic 
component. 

We note that the Chen et al. power-law gives J v = 2.7 x 10 -25 cgs at the Hell edge 
(0.054 keV). This is only slightly smaller than the theoretical Haardt & Madau value of 
3.2 x 10~ 25 cgs at 4 ryd. We therefore assume that the Chen et al. power-law extends down 
to 4 ryd. For lower energies we adopt a second power-law that matches Chen et al. at 4 ryd, 
and Haardt & Madau at the Lyman limit. This gives, 



for vjvQ > 4, where J vo = 2.0 x 10~ 23 cgs and u is the Lyman limit frequency. Most of the 
ionizing energy is emitted near the Lyman limit, although it is evident that the total energy 
in the Chen et al. field (with vJ v ~ z/ - 54 ) increases with photon energy (see Fig. 5). The 
Chen et al. field is truncated at a turn-over energy near 100 keV (see Hasinger 2000). 

Assuming the Vogel et al. (1995) Ha limit, and assuming that the ionizing field varies 
with frequency as given by equations (Bl) and (B2), we derive an upper limit of J UQ < 
7.0 x 10~ 23 cgs for the metagalactic field intensity at the Lyman limit. The upper limit 
(shown in Fig. 8) is 3.5 times larger than the value computed by Haardt and Madau. 




(Bl) 



for 1 < v/vq < 4, and 




(B2) 
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We now consider the FUV and optical extragalactic background. This background is 
generally believed to be produced by the integrated light of faint and distant galaxies. Martin 
and Bowyer (1989) and Martin, Hurwitz, & Bowyer (1991) carried out sounding rocket and 
space-shuttle observations, and detected a 1600 A background flux F* 600 = 40 ± 10 photons 
cm~ 2 s _1 A -1 sr _1 , corresponding to J v = 4.2 ± 1.1 x 10~ 22 cgs at 1600 A. Armand, Milliard 
k Deharveng (1994) found that at 2000 A J v lies in the range of 5.3 x 10~ 22 to 1.7 x 10~ 21 
cgs, consistent with the Martin & Bowyer observations. 

Bernstein (1999) (see also Bernstein, Freedman & Madore 2002) has recently reported 
the first (HST) detections of the extragalactic background light at optical wavelengths. She 
finds that J v equals 1.2 ± 0.57 x 10" 20 , 2.8 ± 0.8 x 10~ 20 , and 4.9 ± 1.3 x 10~ 20 cgs, at 3000, 
5500, and 8000 A respectively. The measured fluxes are consistent with the upper limit of 
1.2 x 10" 19 cgs at 4000 A previously set by Matilla (1989). 

We fit the FUV and optical measurements with a two-piece power-law expression 



for u/uq < 0.3, where again J uo — 2 x 10~ 23 cgs. This fit is consistent with the FUV and 
optical observations within the errors, (see Fig. 5). At 1.25 /xm equation (B4) gives 1.1 x 10~ 19 
cgs, consistent with the (COBE) upper limit of 2.8 x 10~ 19 cgs found by Dwek & Arendt 
(1998) at 1.25 /im. We assume that equation (B4) remains valid down to the H~ threshold 
at 1.65 /im. 

For our fit to the metagalactic field the mean intensity in the 6-13.6 eV FUV range is 
6.7 x 10~ 6 ergs cm -2 s _1 . This is 0.42% of the mean FUV intensity in the solar neighborhood 
as estimated by Habing (1968). 

At what distance from the Galaxy does the metagalactic background dominate the 
radiation field? The Galaxy produces Lyc photons (mainly in OB associations) at a rate of 
Qca\ ~ 2.6 x 10 53 photons s _1 . A small but uncertain fraction, / csc ~ 5%, of this radiation 
may escape the optically thick disk (Dove & Shull 1994; Bland- Hawthorn & Maloney 1999; 
Dove, Shull & Ferrara 2000). For / csc = 5%, the Galactic Lyc flux at the surface of an opaque 
spherical cloud (on the side facing the Galaxy), Qcaifcsc/^d 2 , equals the metagalactic Lyc 
flux of 7rJ* = 3.2 x 10 3 photons cm -2 s™ 1 at a distance d = 180 kpc. The Galaxy may also 
be a source of soft X-rays emitted by cooling hot gas in supernova remnants (Slavin, McKee 




(B3) 



for 1 > v/v Q > 0.3 (i.e. from 912 to 3000 A), and 




(B4) 
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& Hollenbach 2000). These authors estimate that for photon energies greater than 100 eV a 
soft X-ray luminosity of 2.2 x 10 39 erg s _1 escapes the Galaxy. This implies that the Galactic 
and metagalactic X-ray fluxes (between 0.1 and 1 keV) are equal at d — 16 kpc. It is most 
straightforward to estimate the non-ionizing (IR through UV) luminosity of the Galaxy via 
direct observations of an analogous galaxy. For example, Kennicutt (2001) has noted that 
the SBbc galaxy NGC 3992 is very similar to the Galaxy in size, shape, luminosity, and star- 
formation properties. The distance to NGC 3992 is 15 Mpc, and its IR (1.65/xm), optical 
(5500 A), and FUV (1600 A) luminosities are equal to 2.2 x 10 29 , 1.5 x 10 29 and 2.0 x 10 27 erg 
s _1 Hz -1 respectively (data from the NED database 17 ). Assuming these values for the Milky 
Way, the non-ionizing Galactic and metagalactic fluxes are equal at ~ 100 kpc. Given the 
above estimates we conclude that the metagalactic field begins to dominate the ionization 
and heating of the clouds at a distance of 100 to 200 kpc from the Galaxy. 



1 7 nedwww. ipac . caltech.edu 
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phase diagrams for Nm between 10 18 and 10 21 cm -2 , (b) the critical pressures 
.Pmin and P max as functions of the primary ionization rate ( p (s -1 ). The upper 
axis in panel (b) indicates the associated HI shielding columns (cm -2 ) for a 
plane-parallel geometry 71 
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6 M WM vs. M vir diagrams for 10 4 K WM clouds in median (a = 0) Burkert 
halos, and a bounding pressure -Prim — 1 cm" 3 K: (a) total cloud radius r-w/H, 
(b) the ratio, Mwm/^dm, of the WM cloud mass and enclosed DM mass, (c) 
WNM radius %nm, (d) total HI gas mass, (e) peak HI column density, (f) 
central pressure. The dotted curves in panels (c) and (e) indicate the locus 
of critical clouds. The thick solid curves in panels (b), (c) and (f) indicate 
marginally bound clouds. The dashed curves in panel (f ) indicate the "multi- 
phased strip" within which multi-phased behavior is possible in the cloud cores. 72 

7 Mwm vs. M vir diagrams for median (a = 0) Burkert halos, and Phim = 10 
cm" 3 K 73 

8 M WM vs. M vir diagrams for median (a = 0) Burkert halos, and Phim = 100 
cm" 3 K 74 

9 M WM vs. M vir diagrams for —Aa NFW halos, and Phim = 1 cm" 3 K 75 

10 M WM vs. M vir diagrams for — 4cr NFW halos, and a bounding pressure Phim = 

10 cm" 3 K 76 

11 M WM vs. M vir diagrams for -4cr NFW halos, and P H im = 100 cm" 3 K. . . . 77 

12 Median M vir = 2 x 10 9 M (or v s = 30 km s" 1 ) Burkert halo model for the 
dwarf galaxies Leo A and Sag DIG: (a) WM and HI gas densities (solid curves), 
and the sum of the DM and stellar densities (dashed curve) as functions of 
the cloud radius r, (b) projected HI column density TVhi, and the average HI 
column, < N m >, within r, (c) enclosed WM and HI gas masses (solid), and 
enclosed DM plus stellar masses (dashed), (d) Ha recombination line surface 
brightness, and (e) WM gas pressure, Pwm, and the critical pressures P min 

and P max , and (f) total (solid) and primary (dashed) hydrogen ionization rates. 78 

13 Circumgalactic CHVC model. Median Burkert halo, with M vir = 1 x 1O 8 M 
(or v s = 12 km s" 1 ), and Phim A = 50 cm" 3 K. (a) WM and HI gas densities 
(solid curves), and the DM density (dashed curve) as functions of the cloud 
radius r, (b) projected HI column density iVm, and the average HI column 
within r, the dotted curve is the HI column density for an exponential HI gas 
density distribution (see text), (c) enclosed WM and HI gas masses (solid), 
and DM mass (dashed), (d) Ha recombination line surface brightness, (e) 
WM gas pressure, Pwm, and the critical pressures P m j n and P max , and (f) 
total (solid) and primary (dashed) hydrogen ionization rates 79 

14 Galactic HIM pressure for a T H im = 2 x 10 6 K corona 80 
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15 Extragalactic CHVC model. -4<7 Burkert halo, with M vir = 1 x 1O 8 M (or 
v s = 12 km s -1 ), and Pmm/k = 20 cm -3 K. (a) WM and HI gas density dis- 
tributions (solid curves), and the DM densities (dashed curve), (b) projected 
HI column density Nm, and the average HI column within r, the dotted curve 
is the HI column density for an exponential HI gas density distribution (see 
text), (c) enclosed WM and HI gas masses (solid), and enclosed DM plus stel- 
lar masses (dashed), (d) Ha recombination line surface brightness, (e) WM 
gas pressure, Pwm, and the critical pressures P min and P max , and (f) total 
(solid) and primary (dashed) hydrogen ionization rates 81 

16 The small-rc 1/e gas scale height r gas for a 10 4 K WNM in NFW halos: (a) 
r gas for median, +3(7, and — 3<r halos as a function of v s , (b) r gas as a function 
of M vir , (c) level curves for r gas as a function of a and v s , and (d) level curves 
for r gas as a function of a and M vir . The dashed curves in panels (a) and (b) 
display the scale height as given by the analytic expressions in Table 6. . . . 82 

17 The small- x 1/e gas scale height r gas for a 10 4 K WNM in Burkert halos: (a) 
r gas for median, +3<r, and — 3a halos as a function of v s , (b) r gas as a function 
of M vir , (c) level curves for r gas as a function of a and v s , and (d) level curves 
for r gas as a function of a and M vir . The dashed curves in panels (a) and (b) 
display the scale height as given by the analytic expressions in Table 6. . . . 83 
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Table 1. Dwarf Galaxies: Observations 





unit 


Leo A 


Sag DIG 


distance 


kpc 


690 


1061 


optical core radius 


kpc 


0.185 


0.125 


L v 


10 7 L Q 


3.0 


6.9 


single dish 21 cm flux 


Jy km s _1 


69 


33 


VLA 21 cm flux 


Jy km s _1 


65 


30 


WNM velocity dispersion 


km s _1 


9.0 


9.5 


WNM mass fraction 




0.80 


0.83 


AWak (WNM) 


cm -2 


2.2 x 10 21 


1.0 x 10 21 


# g as (projected) 


arcsec 


160 


70 


r N m 


kpc 


0.5 


0.4 


^Hlmin 


cm -2 


2 x 10 19 


5 x 10 18 


r HImin 


kpc 


1.1 


1.6 
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Table 2. Dwarf Model: median Burkert halo 



Halo Parameters 



Xvir 




25.6 


M vir 


M 


2.0 x 10 9 


'"vir 


kpc 


32.5 


n ds 


amu cm' 3 


1.24 


r s 


kpc 


1.27 


v s 


km s _1 


30.0 


M ds 


M 


2.67 x 10 s 


r gas (small-x) 


kpc 


0.49 


Gas Input Parameters 


Pmu/k 


cm 3 K 


1.0 


M W M 


M 


2.2 x 10 7 


M WM /M vir 




1.1 x 10~ 2 


Gas Output Parameters 


?"W/H 


kpc 


4.44 


r gas (numerical) 


kpc 


0.56 




kpc 


0.59 


^WNM 


kpc 


1.39 


^Hlpcak 


cm -2 


1.42 x 10 21 


M m 


M 


1.34 x 10 7 


W/T 




4.4 


Pol PmM 




4.97 x 10 3 


™H,0 


cm -3 


0.46 



Table 3. CHVC Circumgalactic median Burkert halo 



Halo Parameters 







32.5 


M vir 


M 


1.0 x 10 s 


'"vir 


kpc 


12.0 


nds 


amu cm' 3 


2.33 


r s 


kpc 


0.37 


v s 


km s" 1 


11.9 


M ds 


M 


1.22 x 10 7 


r gas (small-x) 


kpc 


0.36 


Gas Input Parameters 


PmM/k 


cm 3 K 


50.0 


M W m 


M 


1.1 x 10 6 


M WM /M vir 




1.1 x 10" 2 


Gas Output Parameters 


rw/H 


kpc 


1.30 


r gas (numerical) 


kpc 


0.55 


r N m 


kpc 


0.50 


?~WNM 


kpc 


0.67 


r CNM 


kpc 


0.33 


Nmpeak 


cm" 2 


5.10 x 10 19 


M m 


M 


2.81 x 10 5 


W/T 




1.0 


Po/ -PfflM 




4.7 




cm" 3 


2.08 x 10" 2 
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Table 4. CHVC extragalactic Model: -4a Burkert halo 



Halo Parameters 







8.1 


M vir 


M 


3.2 x 10 s 


'"vir 


kpc 


17.7 


nds 


amu cm' 3 


7.0 x 10" 2 


r s 


kpc 


2.17 


Vs 


km s" 1 


12.1 


M ds 


M 


7.43 x 10 7 


r gas (small-x) 


kpc 


2.1 


Gas Input Parameters 


PmM/k 


cm 3 K 


20.0 


M W m 


M 


7.1 x 10 7 


M WM /M vir 




0.22 


Gas Output Parameters 


rw/H 


kpc 


7.15 


r gas (numerical) 


kpc 


3.07 




kpc 


2.63 


*"WNM 


kpc 


3.40 


r WNM/CNM 


kpc 


1.16 


^Vfflpcak 


cm" 2 


9.99 x 10 19 


M H I 


M 


1.44 x 10 7 


w/r 




1.0 


Po/ -PfflM 




4.4 


™H,0 


cm" 3 


7.65 x 10" 3 
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Table 5: Halo Profiles 



NFW 



x(l+x) 



/ M = 3[ln(l+x)- T ^] 



ln(l+a:) 

x 
\2 



/gas = e" 3 ^^ (l + x)^^ 



r/x 



Burkert 



(l+x)(l+x 2 ) 



f — 3 

JM ~ 2 


rm(i+* 2 ) + ln(1 + x) tan -i J 






f<P = 1' 


'tan- 1 x)(l + i) - [ln(l + x)](l + ±) + |[ln(l 










(3/2)K/c 9 )2 


/gas = 


e -(l+l/a;)tan- 1 a :^ 1 + ^ (1+1/x) ( J + ^ (l/2)(l/x-l) 



Tabic 6: Small-i Formulae 



NFW 



/ gas = exp (-l^fz 

X erne 



- 2fi 

3 u? 



' gas 3 v 2 ' s 2TrGm H n d3 r a 

r gas = 0.19 x lO- - 186 ^^ - 681 kpc 

r gas = 0.17 x 3 kpc 

Wwm = 47m H , ra H ?- gas Jq 00 y 2 e- y dy = STrn^omnr^ 

JV WM (6) = 2n H , r ga J co e-V( 62 +» 2 )d 2 y = 2n H , r gas 6^i(6) 

^WM.peak = JV(0) = 2n H ,oVs 

EM = 2n 2 H0 r gBS / °° e- 2 »cfy = n H , r gas 

n H ,ocri t = 9.06 x 10- 3 x 10 ' 093 " J 3 * 1/2 c^!? 6 341 cm- 3 

n H ,oc rit = 9.58 x 10- 3 x lO - 108 ^* 172 ^^ 198 cm- 3 



Burkert 

/gas = 6xp ^— 2 ^ 

•kgas — V ^ „ a — Ts 

r r a _ [ 3=1 1 1/2 
'"gas = V2c 9 ^- = ^ 2T Gm H n da J 

r g as = 0.41 x 10-°- 18fe c 96 ^° 6 319 kpc 



Mi 



as = 0.44 x 10-°- 17s - Cfl6 M v i™ kpc 
r wM = 47TO H ,om H r 3 as J °° y 2 e~ y2 dy = 7r 3/2 n H , 77i H ?- gas 
jVwm(6) = 2n H ,oV / °° e^+^dy = 2n H ,ove^ i 

jVwM.peak = JV(O) = v^^H.oVs 

™ = 2«h,oVs Jo 00 e _%2 <iy = vf n|,oV» 
n H ,oc rit = 5.51 x 10- 3 x 10°- M3 -J3 1/2 Cfl6 1/2 ^- 6 - 186 cm" 
n H ,oc rit = 5.32- 3 x 10°- 088 -J* 1/2 Cfl6 1/2 



-0.050 



cm 



Note. — In these expressions, v s , r s , and n^s are the halo scale parameters (see §2.1), and c s is the gas 
sound speed. v s e = v s /10 km s _1 , c. g n = c 9 /10 km s _1 , M Tir 8 = M vir /10 8 M , and J| = J*/10 3 photons 
cm~ 2 s _1 sr _1 . ifi is the first order modified Bessel function. 



